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Abstract 

Starting from the global parametrized post-Newtonian (PPN) reference sys- 
tem with two PPN parameters 7 and (3 we consider a space-bounded sub- 



^ ■ system of matter and construct a local reference system for that subsystem 

<^ . 

fN . in which the influence of external masses reduces to tidal effects. Both the 

(<— ^ ■ metric tensor of the local PPN reference system in the first post-Newtonian 

0^ ■ approximation as well as the coordinate transformations between the global 

O^i PPN reference system and the local one are constructed in explicit form. The 

bJ!)' terms proportional to ?? = 4/3 — 7 — 3 reflecting a violation of the equivalence 



principle are discussed in detail. 



j^ . We suggest an empirical definition of multipole moments which are in- 

tended to play the same role in PPN celestial mechanics as the Blanchet- 
Damour moments in General Relativity. We also show that the tidal grav- 
itational field as seen in the local PPN reference system can be expanded 
into powers of local coordinates similarly to the tidal expansion in General 
Relativity. 

Starting with the metric tensor in the local PPN reference system we 
derive translational equations of motion of a test particle (an Earth satellite) 
in that system. The translational and rotational equations of motion for center 
of mass and spin of each of A^ extended massive bodies possessing arbitrary 

1 



multipole structure are derived. All equations of motion are presented also in 

the form of multipole expansions. Several interesting features of the equations 

are discussed. 

As an application of the general equations of motion a monopole-spin 

dipole model is considered and the known PPN equations of motion of mass 

monopoles with spins are rederived. For the first time, these equations are 

derived in a self-consistent manner which does not require any additional 

assumptions about the behavior of bodies such as secular stationarity. 
PACS numbers: 04.25.Nx, 04.80. Cc, 95.10.Ce 



Typeset using REVT^ 



Contents 



[rr 



III 



[V 



introduction 



"Notations 



Metric tensors and coordinate transformations 



A The PPN metric tensor in the global reference system 



18 

20 

20 



B Energy-momentum tensor for an ideal fiuid| 23 

25 



C The PPN metric tensor in the local geocentric system 



D Coordinate transformations between the global and local reference systems] 27 



IV Matching of the global and local PPJN metric tensors] 27 

A Matching of order C>(c~^) in qoo\ 28 

B Matching of order 0(c~^) in gjj 29 

30 



C Matching of order 0{c '^) in~^ 



D Matching of order 0{c "^) in gpol 32 

E Energy- momentum tensor in the local reference system 



Other forms of the local metric 



VI Multipole expansions of the internal potentials 



Vli Tidal expansions of the external potentials 



36 



37 



39 



41 



Vlli Ttanslational equations of motion of a test particle in the local reference 
systemi 



44 



A The equations of motion in closed form] 44 



B Multipole expansions of the equations of motion] 45 



iX lYanslational and rotational equations of motion of J\ extended bodie^ 47 



A Local equations of motion: T t)^ = 0] 



48 



X 



7179- 



B Local equations of motion: T /^ = 



C Local equations of motion: Egbc X T n = 



D Multipole expansions of the relations between the Tolman and BD masses 



and mass dipole^ 



E Multipole-expanded local translational equations of motion 



F Multipole-expanded rotational equations of motion 



G Mass monopole-spin dipole approximation 



Outlook and conclusions 



Acknowledgment^ 



References 



50 
51 

52 
53 
55 
56 

62 

65 

67 



I. INTRODUCTION 

Einstein's theory of gravity (the so-called General Theory of Relativity, GRT) in the last 
several decades has evolved far from being predominantly a playground for mathematicians. 
Not only that the field of experimental gravity has expanded and the number of experimental 
tests of GRT increased drastically, but certain results from GRT are taken for granted and 
implemented in software for solving technologically oriented problems. A good example 
is Very Long Baseline Interferometry (VLBI) that determines baselines on the surface of 
the Earth, and the position of the rotation pole with subcentimeter accuracy, the length 
of the day at the fraction of a millisecond level, and relative angles between two remote 
radio sources with a precision better than a milliarcsecond. In the VLBI analysis software 
relativistic effects like the gravitational time delay are routinely taken into account. Note 
that the gravitational light deflection amounts to 1.75" at the limb of the Sun and it decreases 
with the inverse impact parameter of the (unperturbed) light-ray from the Sun. This implies 
that the light deflection (more precisely the gravitational time delay) in the field of the Sun 
is measurable even for radio sources lying closely to 180° from the Sun |l|J^. 

For the analysis of high precision optical astrometric measurements at the milliarcsecond 
level (e.g., HIPPARCOS) GRT became an indispensable tool [^. For planned space inter- 
ferometry with microarcsecond accuracies (DIVA, GAIA, SIM, FAME, etc.) relativity will 
play a crucial role and the whole process of observation (essentially, following the photons 
from the emission region onto the CCD-sensor) has to be formulated within the framework 
of Einstein's theory of gravity 0,^. 

Astrometry is one important field of application of GRT, relativistic astrophysics (e.g., 
the physics of white dwarf stars, neutron stars, black holes or quasars, gravitational wave 
physics, cosmology) and celestial mechanics are further ones. Techniques used to gather 
information about the dynamics of astronomical bodies comprise the timing of millisecond 
pulsars at the sub-microsecond level, laser ranging to selected artificial satellites such as 
LAGEOS and to retro corner reflectors that have been placed on the lunar surface at the 



centimeter level and radar ranging to planets or spacecrafts (with, e.g., a few meters accuracy 
for the Viking landers on Mars). 

Meanwhile it is well recognized that GRT enters every astronomical observation or mea- 
suring technique at a certain level of accuracy. Such a critical level has been reached for a 
variety of geodetic measurements including GPS and a long time ago already for the proce- 
dure to establish a terrestrial time scale (TT, TAI). In 1991 the International Astronomical 
Union and the International Union for Geodesy and Geophysics adopted recommendations 
implying the use of GRT in modeling of modern high-accuracy observations P] . The Inter- 
national Earth Rotation Service responsible for precise monitoring of the Earth orientation 
parameters has also adopted relativistic models for several kinds of astronomical observations 

i- 

For practical applications one faces the problem of solving Einstein's field equations for 
complicated situations. Methods of numerical relativity might be employed in the future 
not only for problems related with gravity wave generation (e.g., for the description of 
coalescing binaries) or cosmology, but also for celestial mechanical problems in the solar 
system. However, the classical formulation of celestial mechanics is simplified considerably 
by the introduction of collective variables like mass-mult ipole moments, center of mass, etc. 
and one clearly wants to keep these advantages in a relativistic formulation. Note that the 
use of such collective variables is related with expansions in terms of a 'geometric coupling 
parameter' a = L/R , measuring the ratio between a characteristic linear dimension, L, of 
the bodies and a characteristic separation, R, between them. 

Already soon after Einstein published his first papers on his theory of gravity, Einstein 
[§|, Droste 0, De Sitter [1^, and Lorentz and Droste [0 devised an approximation method 



(called "post-Newtonian") which allowed them to compare General Relativity with New- 
ton's theory of gravity, and to predict several "relativistic effects" in celestial mechanics, 
such as the relativistic advance of the perihelion of planets, and the relativistic precession 
of the lunar orbit: this post-Newtonian approach to general relativistic celestial mechanics 
was subsequently developed by many authors, notably by Fock [|^, Papapetrou |jl3|, Chan- 
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drasekhar and colleagues [1I4|-[T6[], Caporali [jT^, Grislicliuk and Kopejkin [Q and many 



others (for a review of the development of the problem of motion in General Relativity see, 



e.g., Damour |T9|). Basically the post-Newtonian approach is a slow-motion and weak- field 
approximation to Einstein's theory of gravity. One introduces one or several dimensionless 
parameters for the A^-body problem: /3e = forb/c (orbital velocities). Pi = fjnt/c (internal 
or intrinsic rotational velocities), He = GM/{c^R) (external gravitational potential) and 
Ki = GM /{(?L) (internal gravitational potential) if M is the mass of a body with dimension 
L and i? is a typical distance to other bodies. Usually in the post-Newtonian approximation 
one assumes /3f < /3g ~ Ke < Kj = e ^ 1 . E.g., e < 10~^ everywhere in the solar system. 

All of the applications mentioned above require a correspondingly accurate relativistic 
theory of the gravitational iV-body problem consisting of 

(i) the external problem: to determine the motion of the centers of mass of the A^ bodies 
and 

(ii) the internal problem: to determine the motion of each body around its center of 
mass. 

It is a remarkable fact of Newtonian celestial mechanics that these two sub-problems are 
coupled only very weakly. The first reason is the fact that the ellipticities (dimensionless 
mass-quadrupole moments or J2-values) and higher mass-multipole moments of solar system 
bodies are all very small. Therefore the motions of the center of masses are dominated by 
the (approximately constant) masses of the bodies independent of internal motion, and the 
couplings of higher mass-multipole moments to higher-order derivatives of the external grav- 
itational potential (i.e., couplings to the external curvature tensor in a geometric language) 
are very small. On the other hand the local environment of a body is influenced by external 
bodies only through the action of tidal forces that grow linearly with the distance from the 
body's center of mass and vary with L according to {R/L)^. This remarkable decoupling is 
sometimes expressed by the phrase "effacement of internal degrees of freedom in the global 
problem and of the external world in the local system" . 



A second reason for this effacement is the principle of equivalence. In its weak form, 
applicable to Newtonian celestial mechanics, it asserts the equivalence of inertial mass with 
gravitational mass that has been tested by experiments with a relative precision of about 
10"^^ (see e.g.. Will |I| for more detail). 

As explained in great detail, e.g., in p]|-p2[], the treatments of both sub-problems in the 
classical post-Newtonian approach are unsatisfactory. One usually considers the A^-body 
problem in some asymptotically flat space-time that can be covered with one single global 
coordinate system, x'^ = {ct, x, y, z) = {ct, x*), i = 1, 2, 3. In the classical formulation of post- 
Newtonian theory concepts like "center of mass", "mass-multipole moments" (including the 
mass) or "mass-centered coordinates" are defined in this global coordinate system. E.g., 
for applications in the solar system the global coordinates usually will be barycentric ones 
and geocentric spatial coordinates classically were defined by X* = x* — x|;(t) , where x|;(t) 
denotes the global coordinates of the "geocenter". Such definitions do not define a useful 
mass-centered frame in General Relativity, in the sense that it does not efface the external 
gravitational field down to tidal effects. In such a geocentric coordinate system relativistic 
effects are of order /3g ~ 10~^ and not of order Kj = (GMe/c^Re) ~ 7 • 10~^° as one would 
expect. Clearly, the classical post-Newtonian theory tries to fit GRT into purely Newtonian 
framework and bears traces of Newton's absolute space and time. 

In the recent decade the relativistic theory of astronomical reference systems in the 
framework of General Relativity has been considered in a detailed manner by a number of 
authors [p^ , p^ , pO| , pl| , p^ , p^ , p^ , p^^ -p9| . These investigations were motivated both by theoret- 



ical interest to construct a proper reference system for a massive material subsystem and 
by practical requirements of astronomy [p0| -|3^ dictated by increased accuracy levels. The 



two most advanced and complete approaches to construct a physically adequate local ref- 



erence system in General Relativity are the Brumberg-Kopeikin approach p0|j21|j25|.26.29 



and the Damour-Soffel-Xu one p2|j27| , |28|j33[| . For the gravitational A^-body problem both 
formalisms introduce a total of A^ -|- 1 different coordinate systems: one set of global coordi- 
nates x^ = {ct, X*) and one set of local coordinates X"' = {cT, X") for each body comoving 
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with the body under consideration. The relation between global and local (e.g., geocentric) 
coordinates describing one and the same event is assumed to be of the general form 

x^{X'^) = a;^(T) + e':,{T)X'^ + i^{T, X), (1.1) 

where ^^ is at least quadratic in the local space coordinate X" and x^{T) describes the 
motion of a central worldline of the Earth (later chosen as the center of mass). Equivalently 
this transformation is described by 

X^{xn = X^{t) + EtitYj, + S'^(i, ^e) , (1.2) 

with r^(t) = X* — a;^(t) and S" is assumed to be at least quadratic in r^. Note that in 
passing from one form to the other one has to bear in mind that x^(T) and X^(t) de- 



scribe two different points on the central geocentric worldline p5| , |33| . Suitable effacement 
conditions in the X" system lead to constraints on the functions appearing in the coordi- 
nate transformations. Note, that these transformations involve space and time coordinates; 
they present generalized Lorentz transformations. Only by such transformations apparent 
Lorentz-contraction effects etc. can be avoided in the geocentric coordinate system. 

Principle features and results of the theory of local reference system in the framework of 
General Relativity are: 

- a description of the metric tensor g^^, in the global system with only two gravitational 
potentials {w, w'^) as functions of x^; similarly, a description of Gap in some local system 
with only two potentials {W, W) ||2^ ; 



- a new theory of astronomical reference frames, including the transformation rules for 
coordinates p0|- p2| , p7| , |29| and metric potentials p7|| ; 



- an improved description, with full post-Newtonian accuracy, of the gravitational struc- 
ture of each body by means of a set of multipole moments, (M^, S"^), of a body A (e.g., 
the Earth) which are linked in an operational way to what can be observed in the local 
gravitational environment of A P^J27[ ]; 
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- a description of the influence of the external world in a local frame by means of some 
suitably defined tidal moments p0|- p2| , p7| , p9| ; 

- translational [^,0 and rotational ||2^ equations of motion with full post-Newtonian 
accuracy and inclusion of all multipole moments; 

- physical adequate equations of motion of a test particle (e.g., Earth satellite) in the 
local reference system [pSJ^pOJ ; 

- physically adequate relativistic models for many kinds of observations (VLBl ||3^36 



high-accuracy positional observations |20.1], various methods of remote clock synchroniza- 



tion 37 , etc. 



Despite the fact that GRT so far passed all experimental tests with flying colours (e.g., 
Misner et al. |38], Will [ffl]) people never lost interest in alternative theories of gravity. 



This is one reason why a parametrized post- Newtonian formalism (the PPN formalism) 
was developed to cover the post Newtonian limits of a whole class of metric theories of 
gravity including GRT. A primitive version of such a formalism was devised by Eddington 
P^ and Robertson [^ who introduced the space curvature parameter 7 and the non- 
linearity parameter (3 that are both equal to one (/3 = 7 = 1) in GRT, but may take 
different values in alternative theories of gravity. Later Schiff [^ included terms from the 
rotation of bodies (gravito- magnetic terms). The classical PPN-formalism was developed 
mainly by Kenneth Nordtvedt in the late 60s and early 70s in collaboration with Clifford 
Will (Nordtvedt [H,!!, Nordtvedt and Will H, Will and Nordtvedt H). More PPN- 
parameters have been added to cover the post-Newtonian limits of a variety of scalar-tensor 
theories like the Brans-Dicke theory, vector tensor theories, bimetric theories, etc. The 
choice of PPN-parameters has changed over the years until some standard was achieved 
that can be found in Will |l[]. Here, a total of ten PPN-parameters is introduced not only 
to distinguish between the various post-Newtonian limits of different theories of gravity, but 
also to provide a tool to characterize the precision of relevant gravitational experiments. In 
the following we will restrict ourselves to the Eddington- Robertson parameters (3 and 7. The 
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other PPN-parameters that are related with possible preferred locations, preferred frames 
or a violation of conservation of total momentum will not be considered here. This is not 
a crucial restriction for our approach that can be extended in a straightforward manner to 
any choice of PPN parameters. However, considering that the two parameters are the most 
important ones we confine ourselves to (3 and 7 only. 

This Eddingt on- Robert son- Schiff Nordtvedt-Will PPN-formalism is formulated in the 
spirit of classical post-Newtonian framework described above. Nevertheless, it was highly 
successful and lead to new predictions, e.g., related with a breakdown of the strong equiv- 
alence principle. Let us summarize briefly the most interesting results from the classical 
PPN-formalism that one would expect to be valid even in an improved version. 

1. In every metric theory of gravity a test body of negligible extension and negligible 
gravitational self-energy will follow a geodesic in a suitably defined external metric g^^^. 
From the usual form of the PPN metric one finds that the geodetic acceleration of a body 
E is given by 

r* 1 

«£,geodesic = - ^ G Mb — h — (. . .) (1.3) 

with suitably chosen masses Mb of gravitating bodies and r^^ = x|; — x^. Note that this 
acceleration (or the geodetic equation) follows from an action integral for the worldline of 
some central point of body E of the form 



1/2 



Cdesic = -Mec j [-g'^^dx^dx''\ . (1.4) 

The corresponding A^-body Lagrangian takes the form [^ (r^^ = x^ — x^) 

^geodesic- I^MaC j^i ^^^ g ^4 j 



^ 1^ ^ GMaMb 

2 A B^A ^AB 



1 + (27 + 1)^-(2/?-1)E|'^' 



^-(2/?-l)E^F^ 



11 1 

-^(47 + 3) va ■ vb - -- (v^ ■ Tab) (vb ■ tab) ^- 



2c2 2c2 , ^^ 



:i.5) 



2. In general a body E with non negligible gravitational self-energy VLe will not follow a 



geodesic because of a violation of the Strong Equivalence Principle (e.g., Nordtvedt |^ , ^ 
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Will |]l|). In that case the value for the gravitational mass of E, M^ , will differ from that 
for the inertial mass Me by 

mP =Me + v^e/c^. (1.6) 

Here, r] is the Nordtvedt-parameter. A classical result gives its relation to the Eddington- 
Robertson parameters 

7^ = 4/3-7-3. (1.7) 

This is the well-known Nordtvedt effect leading to a polarization of the lunar orbit around 
the Earth in the direction of the Sun if r^ 7^ 0. The ratio of gravitational to inertial mass has 
an interesting and valuable interpretation. A violation of the Strong Equivalence Principle 
implies that in the local (e.g. geocentric) system the effects from the outside world cannot 
be effaced down to the usual tidal forces. Instead the "effective" gravitational constant, 
Giocai, as experienced in the local _E-frame becomes effectively space dependent 

Giocai = G ("l - r^^^j . (1.8) 

Here, U is the external gravitational potential produced by the outside world, i.e. by all 
gravitating bodies different from E. Since second and higher order derivatives of U give rise 
to the usual tidal forces the effacement breaks down for the zeroth and first derivative of U, 
taken at the geocenter x^;. Considering the gravitational action between two points x and 
x' one has effectively p6[ 



G^local — G 



UJ^e) V 
c2 2 



uyj<.E I I . n 



;i.9) 



where we have replaced Vf/(x£;) by the acceleration a^;. Deriving the acceleration, e.g., of 
a satellite in the vicinity of the Earth by classical means, but with this expression for Giocai 
leads to the above ratio of gravitational to inertial mass of the Earth. 

3. Already in Newtonian celestial mechanics extended bodies in general do not experience 
just the 1/r^ acceleration when interacting with another body at distance r. Instead higher- 
multipole mass moments (C/m, Sim or equivalently Ml in the Cartesian language with / > 2) 
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lead to a deviation of the 'geodetic' free-fall behavior. To lowest order this force of geodetic- 
deviation is given by 

«£,quadrupole = JJ^ ^jfc ^ijU U{^e) = J^ ifk dijk U{^e)- (1-10) 

Here, 

Ifu = j^P*r'Er%<fx (1.11) 

is the Newtonian second-order mass tensor related with the Newtonian trace-free quadrupole 
tensor, M^^, by 

M% = lf,-\lfj,u. (1.12) 



Note that the last equality in ( p.. 10] ) follows from the Laplace equation in the absence of 



gravitational sources. Since second spatial derivatives of the external potential can be inter- 
preted as components of a Newtonian external curvature tensor this quadrupole acceleration 
of body E can be interpreted as resulting from a coupling of the body's second-moment mass 
tensor with the curvature tensor produced by the outside world. For that reason in any met- 
ric theory of gravity one expects the world line of a suitably chosen central point inside body 
E to be influenced by the body's covariant second- moment mass tensor, J^^, via a coupling 
to the external curvature tensor R^^p,^ as described by the action integral (Nordtvedt 



I = -Me c J 



9j^v + irj~Raiipu 1 ^^ dx 



1/2 

(1.13) 



In the Newtonian limit only the spatial components of /^ (as measured in the local geo- 
centric frame) contribute and from the action integral one derives an acceleration due to the 
second-moment mass tensor of 



c^ 1 — 

C^B.quadrupole = ^ ]\/f J^ ^i^jQkO = 'q1\T~ J^ (JijkU \^e) ■ (1-14) 

As in Newton's theory, in GRT one finds that the trace of iR. does not lead to an additional 
acceleration because of the vacuum field equation (i.e., the vanishing of the Ricci-tensor in 
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source-free regions). However, in general for 7 7^ 1 the Ricci-tensor in vacuum does not 
vanish in alternative theories of gravity. In that case there will be an additional acceleration 
due to the extension-curvature coupling proportional to 7 — 1 even for spherical bodies with 

/S = ^/%. (1-15) 

and 

J^= / pWWe^x. (1.16) 

Je 

This implies that the effacement in GRT in the local (geocentric) system does not only 
lead to the ordinary value of the gravitational constant (since r^ = 0), but also prevents 
extension-curvature couplings for spherical bodies (since 7 = 1). 

4. Finally problems related with the spin (intrinsic angular momentum) have been 
treated extensively in the literature. There are two aspects that have to be tackled: (i) 
the influence of the spin upon the translational equation of motion and (ii) the motion of 
the spin vector itself (the rotational motion). Progress for solving the first problem has 
been achieved by Mathisson |^, Papapetrou |]rB[, Corinaldesi and Papapetrou |^ without 



definitive conclusive results for GRT. However, in the first post-Newtonian approximation 
the situation is different. The spin-orbit and spin-spin coupling terms have been derived by 
Brumberg [Q who also discusses the physical implications in detail (see also Kalitzin |5l| 
and Michalska [^). They have also been derived using a one-graviton exchange theory by 
Barker et al. ||5^ and by Borner et al. |^^ with parameters (3 and 7 (see also Barker et al, 
|55|-p8|). Tulczyjew [^ has derived formally the same results using distributional sources. 



As first discussed by Damour [Q , if we require invariance under Lorentz boosts, the spin- 
orbit Lanrangian depends not only on positions and velocities but also on accelerations. 
Within the Lorentz-invariant Lagrangian formalism this is the manifestation of the well- 



known Thomas precession pT|. The Lagrangian for the spin-orbit and spin-spin coupling 



terms for body E takes the form pT|J53| fl] 



14 









3 (Se ■ rijij) (Sb ■ veb) 



>_B ■ 'JB 



:i.i7) 



' EB ' EB 

Already to Newtonian order the motion of the spin is governed by extension-curvature 
couphngs (precession of the equinoxes; e.g., Misner et al. [0), were the quadrupole and 
higher order mass-multipole moments couple to the external curvature tensor. In the absence 
of such tidal forces the spin vector is Fermi-Walker transported (e.g., Misner et al. ||5^) along 
the body's central world line. This Fermi- Walker transport has been extensively discussed in 
the literature. One finds that the spin- vector in the comoving frame precesses with respect to 
remote celestial objects (more precisely with respect to the local coordinate induced tetrad, 
see e.g. |],|1), i.e. 

-^ = f2xS (1.18) 

dr 

with 

IvxQ IVxW / l\vxVt/ ,^ ^^, 

Here, v is the global coordinate velocity of the rotating body (gyroscope) , Q are the spatial 
components of the body's four-acceleration and W the three space-time components of the 
metric tensor in the comoving system. The three terms on the right hand side describe 
the well known Thomas precession (first term), Schiff or Lense-Thirring precession (second 
term) and the geodetic or de Sitter- Fokker precession (third term). The extension-curvature 



coupling in the post-Newtonian framework of GRT have been discussed in detail in p8[. For 
recent articles on this coupling in GRT see e.g., Apostolatos ||6^j63| . 



It is the purpose of this article to extend and improve the classical PPN framework with 
parameters (3 and 7 so to cover also the principal results from the Brumberg-Kopeikin and 
DSX formalisms. 
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Both Brumberg-Kopeikin and DSX approaches use the field equations of GRT in order 
to construct the local reference system. However, in the PPN formalism we do not have any 
kind of "parametrized post-Newtonian field equations", but only the metric tensor in the 
global reference system. To construct a physically adequate local reference system means to 
find a suitable coordinate transformation which maps the global barycentric metric tensor 
into a "good" local one having some desired properties. These desired properties in General 
Relativity read 

A. The gravitational field of external bodies is represented only in the form of a relativistic 
tidal potential which is at least of second order in the local spatial coordinates and 
coincides with the usual Newtonian tidal potential in the Newtonian limit; 

B. The internal gravitational field of the subsystem coincides with the gravitational field 
of a corresponding isolated source provided that the tidal infiuence of the external 
matter is neglected. 

These two requirements can simultaneously be satisfied in General Relativity as has been 
shown in the framework of the Brumberg-Kopeikin and DSX formalisms. It is clear that this 
fact is closely related with the validity of the Strong Equivalence Principle in GRT. In the 
PPN formalism due to a possible violation of the Strong Equivalence Principle one expects 
that either property A or B can be satisfied, but not both of them simultaneously. 

Both Brumberg-Kopeikin and DSX formalisms rely essentially upon the field equations of 
general relativity. In contrast to this here we do not want to confine ourselves to a particular 
theory of gravity and use instead the phenomenological framework of PPN formalism. It is 
clear that considering a particular theory of gravity might provide a deeper insight into the 
physical origin of various effects in that particular theory. However, this physical meaning 
would be valid only in that theory of gravity and cannot be "generalized" onto other theories 
of gravity even if they give formally the same post-Newtonian metric. 

The global barycentric PPN metric is assumed to be a solution of the field equations of 
any metric gravity theory covered by the PPN formalism. The construction of a physically 
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adequate ("good") local reference system involves a suitable coordinate transformation which 
converts the global barycentric coordinates into "good" local ones. It is clear that if the 
original barycentric metric represents a solution of the field equations, then the resulting 
geocentric metric is also a solution of the same equations. The metric tensor of the local 
reference system is sufficient to discuss both translational and rotational equations of motion 
in that local reference system since they follow either from the geodetic equations (for test 
bodies and photons) or from the local equations of motion (for extended bodies). Both kinds 
of equations require only the metric tensor for their detailed evaluation. No other possible 
scalar, vector, etc. fields appearing in a particular non-Einsteinian theory of gravity play a 
role in discussing equations of motion in the PPN formalism. This allows us to construct the 
local reference system from some geometrical considerations. An attempt to achieve this has 
been already undertaken by Shahid-Saless and Ashby [Q , who introduced a kind of Fermi 
normal coordinates in a simplified model case of two spherically symmetric nonrotating 
gravitating bodies. However it is clear that Fermi coordinates in the PPN framework lead 
to the same principal difficulties when considering the local subsystem with a complicated 
multipole structure, to which they lead in GRT (see, e.g., P^J2^ ). We use the much more 
powerful and elegant approach of the Brumberg-Kopeikin and DSX formalisms. 

A practical application of the theory of relativistic astronomical reference systems lies in 
the construction physically adequate models for modern high-precision astronomical obser- 
vations (depending on the kind of observations the models involve the coordinate transforma- 
tions between the reference systems, equations of light propagation in both reference systems 
as well as equations of translational and rotational motion of extended bodies). However, on 
the other hand, these observations yield a significant amount of experimental data that can 
be used for testing metric theories of gravity such as General Relativity. The Parametrized 
Post-Newtonian (PPN) formalism |I[ presently can be considered to represent the standard 
framework for such tests in the first post-Newtonian approximation. Therefore, it is quite 
important to have a consistent theory of relativistic astronomical reference systems in the 
PPN formalism. Besides a purely theoretical interest to construct a "good" proper reference 
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system for a massive subsystem (a massive body) in the framework of the PPN formahsm, 
a PPN theory of astronomical reference systems is also important from a rather practical 
point of view, since presently there is a contradiction between the recommendations of the 
International Astronomical Union concerning the relativistic reference system to be used 
for data reduction which are valid only in GRT and the common practice to estimate the 
PPN parameters (at least 7 and (3) from observations. Moreover, the International Earth 
Rotation Service being an official service of the lAU recommends in its lERS Conventions 
^ to estimate the PPN parameters 7 and f3 in routine data processing. 

This paper focuses mainly on problems of reference systems and celestial mechanics. 
More specifically, we present 

- a new PPN theory of astronomical reference systems (Sections |ITl|- 0); 

- a new description of the local gravitational environment of a body with new PPN mass- 
and spin-multipole moments generalizing the Blanchet-Damour moments from |^2[ (Section 
[VID and of the tidal forces (Section |V11| ); 

- new PPN satellite equations of motion (Section |V111[ ). 



- new translational and rotational global equations of motion of A^ extended, rotating 
and gravitationally interacting bodies of arbitrary shape and composition with all multipole 
moments of gravitating bodies with parameters f3 and 7 (Section pC]) ; 

A overall outlook and a summary of our main results are given in Section ^. As expected all 
results that have been listed under 1-4 above are recovered, but now with the help of a new 
and improved PPN-framework. Some partial results have already been published elsewhere 



II. NOTATIONS 

Greek indices a, /?,..., /z, z/, .. . running from to 3 indicate all four space-time com- 
ponents of the corresponding variable. Latin indices a,b, . . . ,i,j, . . . run from 1 to 3 and 



refer to three spatial components of the corresponding variable. We use Einstein's sum- 
mation convention for both types of indices independent of the position of repeated in- 
dices: e.g., X* x* = {x^Y + {x^Y + (x'^)^. We use several special objects: 5^^ = diag(l, 1, 1) 
is the Kronecker delta; eijk is the fully antisymmetic Levi-Civita symbol (£123 = +1); 
rj^i, = diag(— 1, +1, +1, +1) is the Minkowki flat space metric (we use the signature (— +++) 
throughout this paper). Capital latin indices A,B,C, . . . = 1 . . . N (as subscripts or super- 
scripts) refer to the bodies of our A^-body system: e.g., ^e is the value of some variable $ 
related to body E. 

The global coordinates are designated by x^ = (ct,x*), t being coordinate time and 
X* being spatial coordinates of the global reference system. The global metric is gfj_!^{x^) 
and we systematically use lower case latin letters for the quantities belonging to the global 
reference system. As in the DSX papers [0,^,^,Q the space-time indices referring to 



global coordinates are taken from the second part of the greek alphabet /i, i/, . . . and space 
indices are taken from the second part of the latin alphabet i,j, . . .. The local coordinates 
for body E are designated by X° = {cT, X"), T being coordinate time and X"- being spatial 
coordinates of the global reference system. The local metric is GapiX'^) (do not confuse an 
index (3 with the PPN parameter /?). The space-time indices referring to global coordinates 
are taken from the first part of the greek alphabet a, (3, . . . and space indices are taken from 
the first part of the latin alphabet a,b, . . .. 

A comma before an index designates the partial derivative with respect to the corre- 
sponding coordinates: A^^ = dA{t,^)/dx^', A^i = dA{t,^)/dx\ A^^ = M(T,X)/9X", 
A^a = dA{T,X.)/dX°'. For partial derivatives with respect to the coordinate times t and 
T we use A^t = dA(t,:K)/dt, A^t = dA(T,X.)/dT. A semicolon before a space-time index 
denotes the corresponding covariant derivative. Parentheses surrounding a group of indices 
denote symmetrization, e.g., 74(jj) = | {Aij + Aji). Brackets surrounding two indices denote 
antisymmetrization, e.g., Aiijk] = | {Aijk — Aikj). Angle brackets surrounding a group of 
indices or, alternatively, a caret on top of a tensor symbol denote the symmetric trace-free 
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(STF) part of the corresponding object, e.g., Aat = A(ab) = STFatAat = A(afe) - lS''''Acc. 
We shall freely lower and raise spatial indices, e.g. Aa = A"-. For sequences of spatial indices 



we shall use multiindices in the same way as it is done, e.g., in [^: a spatial multi-index 
containing / indices is denoted by L {K for k indices, etc.): L = ii . . . ii ii L refers to global 
coordinates and L = ai ... ai in the case of local ones. We use also L — 1 = ii . . . «/_!, etc. A 
multisummation is understood for repeated multiindices: AlBl = Z]ii...i, ^n.-.j; Bi-^,,,i^. For 
a spatial vector f * we denote v^ = v^^ v^'^ . . . f *' . For an L-order partial derivative we denote 
dL = di,... di, . 

III. METRIC TENSORS AND COORDINATE TRANSFORMATIONS 

A. The PPN metric tensor in the global reference system 

We consider a gravitational A^ body system in an asymptotically flat space-time that can 
be covered by one single (global) coordinate system x^ = (ct,x*) with 

lim g^^ = r/^^, (3.1) 

|x| — -oc 
t — const. 

where (/^i, is the metric tensor in the global coordinate system and r/^i, = 
diag(— 1, +1, +1, +1) is the fiat metric tensor of Minkowski space-time. The metric ten- 
sor Qui, in the global PPN reference system is written as 



c2 ^ ' ' d 



2(1 + 7) 



-5\ 



C 



2 



g,, = 6,, (^1 + -7w;(t,x)j + 0{c-^). (3.2) 

In accordance with the classical PPN framework we will assume the metric potentials w and 
w* to obey the equations 

w,u-^w,tt = -4:71 Ga + 0(0-^), (3.3) 



c 



w]j^ = -ATiGa' + 0{c-'), (3.4) 
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where 

cr = i (^tO° + 7T'=' + ^T°°(37-2/?-1)«7)+0(c-^), (3.5) 

cr' = -T°' + 0(c~'). (3.6) 

c 

Here, T'^'^ are the components of energy-niomentum tensor in the global reference system 
and w in ( |3.5| )is needed only to Newtonian order where it coincides with the Newtonian 
potential. Because of requirement ( |3.1D {w^ = {w,w'')) we have 

lim w''{t,x) =0 (3.7) 

|x|— *oc 
t — const. 

and the solution of ( |3.3| )-( p^ ) can be written in the form 

w^{t, ^) = G f p^d'x' + ^G^ (a^{t, x') |x - xV'x' + 0{c-'). (3.8) 

J |x — x'l zc^ oV J 

The expressions above are valid independently of the particular form of the energy- 
momentum tensor. This is no more than a formal way to specify the global PPN metric 
tensor. If matter is described as an ideal fluid or stressed continuum the metric ( ^.2D -( P?B[ ) 
coincides with the version of the PPN formalism described in |I]] as well as with the version 
discussed in [^, provided that only two parameters 7 and (3 are retained in both versions. 
Effectively we consider only those theories of gravity which produce the metric (p.2|) - (|3.6|) 
in the first post-Newtonian approximation for some form of the energy-momentum tensor. 

Eqs. (|3.3| ) - (|3l6|) can be considered as "effective" or "phenomenological" post-Newtonian 
field equations in the global reference system. These equations reduce to the Einstein field 
equations for 7 = /? = 1. It is, however, clear that these phenomenological field equations 
have nothing to do with the actual field equations of a specific theory of gravity covered by 
the PPN formalism. 

The system under study is supposed to be spatially bounded and isolated which means 
that T^" = sufficiently far from the origin. We suppose also that the matter in concentrated 
in spatially separated blobs (bodies) and T^*^ vanishes outside these blobs. One of such blobs 
is selected to define a local subsystem and the matter inside the blob is referred below to 
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as internal matter. All other matter is called external. Any subsystem is characterized by 
a world tube V which encompasses the internal matter. Since the subsystems are supposed 
to be spatially separated we choose the world tube V in such a way that T'^'^ = on the 
boundary of V. In the following we will concentrate upon one of these bodies that we will 
call "Earth" (indicated by subscript E). 

The potentials w'^ are defined by (|3.8| ) as volume integrals over the whole 3-space. There- 
fore, splitting the area of integration into the volume V of the body, for which we want to 
construct a local reference system, and the remaining part of space, we split w'^ into internal 
potentials (potentials of the body under consideration) and external ones (potentials due to 
the other bodies): 

w{t,X.) = WE(t,x) +w{t,X.), 

w\t, x) = ^(t, x) + w\t, x). (3.9) 

The internal potentials w^ = {we, w^^) are defined as 

w^{t,^) = G f p^d'x' + ^G^ /a^(t,x')|x-x'|rfV. (3.10) 

Jv |x — x'l zc' ot'^ Jv 

We adopt a harmonic-like gauge for the global PPN metric tensor. Precisely speaking, 
the global metric tensor satisfies the usual harmonic gauge {g = det{gfj_^)) 

J:r((-9)"V)=o (3.11) 

in case of General Relativity /? = 7 = 1. This requires 

w,t + wii = 0{c-^), (3.12) 

which is also compatible with the definitions of w and w\ and the Newtonian continuity 
equation 

a,t + a:, = Oic'^). (3.13) 

Since the energy-momentum tensor T'^'^ is supposed to vanish on the boundary of the world 
tube V encompassing the subsystem, Eq. ( |3.13| ) allows one to formulate ( ^.12] ) separately 
for internal and external potentials 
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WE,t + W^^i = 0{C 



(3.14) 



w^t + w\ = 0{c 



(3.15) 



The metric ( |3.2| )-( [3^ ) is equivalent to tiie PPN metric in [Q up to a trivial gauge 
transformation. The transformations between our reference system (t,a;*) and the PPN 
reference system in the standard post-Newtonian gauge (tpN, a^p^) reads 



tpN = t - ^ X,t + 0{c 



X 



PN 



X' 



where x is the superpotential, here defined by 



X 



]-G I a{t, x') |x - x'l dV + C(c-2), 



so that 



X,ii = w + 0{c 



(3.16) 



(3.17) 



(3.18) 



The transformation of the results derived below into the standard PPN gauge represents 
no principal difficulty and can be done in exactly the same way as it has been done in Section 
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B. Energy-momentum tensor for an ideal fluid 



Although, as noted above, our formalism is not restricted to any particular form of the 
energy-momentum tensor, the standard PPN formalism (see, e.g., |Q) is elaborated for the 
case of an ideal fluid and it is interesting to specify our formalism for this simplified case. 
The energy-momentum tensor of an ideal fiuid in the global asymptotically fiat reference 
system (t, x*) reads 



TOO ^ p^2 



r 



OJ 



pcv 



l + -[v^ + Il + 2U 



+ 0{c~ 



+ -pv' + 0{c-''), 



T^ = pv'v^ + 6'^p + 0{c- 



(3.19) 
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Here f* = f*(t,x) is the velocity of matter, n(t,x) is the specific energy density (ratio of 
energy density to rest-mass density), p{t,x.) is the isotropic pressure, p(t, x) is the invari- 
ant density (that is, the density invariant under coordinate transformations from (t, x) to 
any (t, x): p(t,x) = p(t, x)). In many cases it is more convenient to express the energy- 
momentum tensor in terms of p*{t, x), 

1 /I 



p*(t,x)=p(t,x)(-^)V2^ = p 



l + ^i^^v^ + 3^U)+0{c-'''^ 



(3.20) 



For the density p* a Newtonian-hke continuity equation holds to post-Newtonian accuracy 

m^-^ilio-^'l-Oic-^ (3.21) 



In terms of p* the energy-momentum tensor reads 



T^OO 
rpOi 



* 2 

p C 



p cv' 



1 /I 



v^ + U-{3j-2)U 



1 + 



1 /I 



c2 V2 



n - (37 - 2) [/ 



+ 0{c- 



c 



p*v'v^ + 6'^p + 0{c-^). 



(3.22) 



Here and above U is the Newtonian-like gravitational potential 

t/(t,x) = G I 



"•''^dV. 



X — X 



For an ideal fluid the potentials w and w^ can then be expressed as 
w = [/ + — f (7 + -) $1 + (1 - 2/3)$2 + $3 + 37$4 + X,tt 



z.- = G/ ^*(^'">^(5"^^ ^V + 0(c-^) 



$3 = 6' 



|x — X' 



p*(t,xOi^^3^, 



X — X' 



P*(t,x0[/(t,x0 ^3^, 



X — X' 



P*(t,x0n(t,x0 3^, 

I /I 

X — x' 



-dV, 



$, = G / ^^^rfV. 



(3.23) 



X — X' 



(3.24) 
(3.25) 
(3.26) 
(3.27) 
(3.28) 
(3.29) 
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Since all the potentials ( |3.23| ), (|3.26|) - (|3.29|) and (|3.17|) are again some volume integrals, 
we can split them into local parts where the integration is performed over a world tube V 
covering the subsystem for which the local RS is being constructed and external ones (the 
integration is performed only outside of V) 

U = Ue + U, 

$, = $f + $,, ^ = 1,2,3,4, 
X = X'' + X- (3.30) 

For example, the definition of Ue reads 

f/E(t,x) = G f ^l^dV. (3.31) 

Jv |x — x'l 

C. The PPN metric tensor in the local geocentric system 

We assume that the metric tensor G^y{T, X) in the local reference system X^ = {cT, X°-) 
has the same functional form as the metric in the global reference system 

Goo = -1 + ^ W^(T,X) - 4/5lV'(T,X) + 0(0"^), 



Goa = -^^^^ W^iT, X) + 0{c 



2(l + 7)^a^y^X) + 0(c-^) 
2 



C3 



Gab = 5ab[l + ^7W^(r,X) ) + 0(c-^), (3.32) 



and the potentials W and W^ entering into ( p.32 ) admit a representation in the form 



WiT, X) = W^iT, X) + QaiT) X" + Wt(T, X) + ^ *(T, X), (3.33) 

W^"(T,X) = H^^(r,X) + ^£,bca(T)X^ + W^(T,X). (3.34) 

Here Wt and W^ are external potentials representing tidal fields of the other bodies of the 
system and are assumed to be (9(X^). The two arbitrary functions Qa and Ca have a clear 
physical meaning which will be discussed below. The local internal gravitational potentials 
are represented by three other functions VFe, W^e ^^"^ ^ ■ The potentials We and W^ are 
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supposed to have the same functional form as their counterparts we and Wg, but expressed 
in local coordinates: 

We = G f S(T, X') ^ d'X' 
Jv |X — X'l 

W^ = G f J:%T,X') \ d'X' + 0{c-'), (3.36) 

jy |X — X I 

S = 1 ("t"" + 7T"" + 4^°° (37 - 2/3 - 1) ly) + 0(c-'), (3.37) 

2^^ = -r°'^ + C(c-2), (3.38) 

c 

where T^'^ is the energy- momentum tensor in the local reference system (T, X"). Finally, the 
function \1'(T, X) represents any possible deviation of the actual internal gravitational field 
as seen in the local reference system from the form (|3.35|) - (|3.36|) . Clearly the appearance of 
\1/ is related with a violation of the Strong Equivalence Principle which makes it impossible 
to satisfy simultaneously properties A and B formulated in Section |I[ By assuming that Wt 
and W| are ~ (9(X^) we assume that property A is satisfied. For that reason property B 
will be violated resulting in the appearance of the additional quantity ^. 

The local metric is also supposed to satisfy the harmonic gauge {Q = det^Gap)) 

_|_ ((_^)1/2g'«;3) ^ (3.39) 

in case of General Relativity (/3 = 7 = 1). This requires 

W,T + W^^ = 0{c-'^). (3.40) 

From (|3.35| ) and (|3.36|) , the Newtonian continuity equation 

Et + S% = 0{c'') (3.41) 

and from the assumption that the energy-momentum tensor (and, therefore, S and E") 
vanishes on the boundary of the world tube V it follows 

26 



W^,T + Wl, = 0{c-^). (3.42) 

From d^ig) and ( |02[ ) one finds 

QaX^ + PFt.t + W^T.a = (^(c"')- (3.43) 

Let us note that at this point there is no way to prove that a local reference system 
with all these properties exists. Unlike the case of General Relativity we do not have field 
equations that would allow us to prove some properties of the local reference system. We 
just formulate the desired properties of the local reference system and then prove that it is 
indeed possible to construct such a system. 

D. Coordinate transformations between the global and local reference systems 

The results of the Brumberg-Kopeikin and DSX formalisms (see, e.g., Theorems 1 and 2 



of [^) allow us to write the transformations between the global and local reference systems 



in the form 



T = t-\[A + 44) + \{b + B% + B'^r\,ri + C(t, x)) + 0(0"^), (3.44) 



c^ ^ — ^ c 

1 /I 



X'' = i?;;. [ri + - [j^vivir'^ + D^'^rl + D^'^' r|rj, j j + 0(c-^), (3.45) 

where r^ = x* — ^^(t), ^^(t) is the coordinates of the origin of the local reference system 
relative to the global one, and f^ = dx^^/dt and a^ = d'^x^^/dt'^ are its velocity and accel- 
eration, respectively. The functions A{t), B{t), B\t), B'^{t), D'^{t), D'^^{t), R^jit) (being 
a rotational (orthogonal) matrix) and C(t,x) ~ 0{r'^) are some unknown functions to be 
determined with the aid of matching of the two metric tensors. 

The fact that the metric ( |3.2| ) can be transformed into the metric ( |3.32| ) with the help of 
the transformations ( |3.44| )-( |3.45|) is not obvious from the beginning and has to be verified. 



IV. MATCHING OF THE GLOBAL AND LOCAL PPN METRIC TENSORS 



Matching of the metric tensors ( p.2|) and ( p.32| ) 
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BX^ f)X^ 

where the Jacobian dX^ jdx^ has to be determined from the coordinate transformation 
( p.44| )- (|3.45| ), yields the exphcit expressions for the unknown functions in the coordinate 
transformations as well as in the local metric. The matching is performed in four steps, 
considering 



A) terms of order 0{c~ ) in ^fQo (e = A = in (|4l|)); 

B) terms of order 0{c~^) in gij {e = i, X = j in ([4.1|)); 

C) terms of order 0{c^^) in goi {e = i, X = in (^)); 

D) terms of order 0{c~'^) in goo (e = X = in (|4l|)). 

At each step of the matching procedure we equate separately four kind of terms in the left- 
and right-hand sides of (|4.1| ): 

1) terms containing the internal potentials of either global or local metric; 

2) external terms independent of local spatial coordinates (that is, the terms which are 
functions of time only) ; 

3) external terms linear with respect to local coordinates X; 

4) external terms of at least second order with respect to local coordinates X. 
This procedure allows us to derive the following results. 

A. Matching of order 0{c^'^) in goo 

1) Internal terms give the Newtonian transformation between the internal potentials in the 
global and local reference systems 

We{T,X) = WE{t,^) + 0{c-^). (4.2) 

28 



It is easy to see that ( [4 .21 ) is in agreement with the definitions of WE(t,x.) and We{u,w) 
(see, (|3iq ) and (p:35|)). 



2) External terms of order \'X.\^ give 

^A=l4 + W(x^). (4.3) 

Here and below for any function $(t, x) we use the shorthand notation $(xe) = ^(t, X£;(t)). 

3) External terms of order |X|^ lead to the Newtonian equations of motion of the center of 
mass of the local system relative to the global reference system 

a^(t) = W,i(xs) - i?^ Qa + 0{c-^). (4.4) 

4) External terms of order |X| , / > 2 yield the expression for the external tidal potential in 
the local metric 

Wt{T, X) = w{t, x) - W(xs) - tZ;,,(xs) r^ + Oic~^). (4.5) 

Here and below for any function $(t, x) the shorthand notation $ j(x£;) denotes d^{t, x)/9x' 
evaluated at x = X£;(t). 

B. Matching of order 0{c^'^) in gij 

1) Internal terms are matched automatically due to ( [4.21) . 

2) External terms of order \X\^ give 

D'^{t) = 5'^-iw{^E)- (4.6) 

3) External terms of order |X|^ give 

D''\t) = ^ 7 {6'^a% + 6'^a^E - ^^'a'^) . (4.7) 
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Let us note here that this expression for D'-'*^ follows from the isotropic form of Qij and Gab 



and its relation to goo and Goo, respectively (this was proven by Theorem 2 of ||22| in the 
case 7 = 1). 

4) External terms of order |X| , / > 2 are matched automatically due to ( [4.5|) . 

C. Matching of order ©(c •^) in goi 

1) Equating internal terms results in the expression for the local vector potential 

W^iT,:S.) = R% {wUt,^)-v'j,WE{t,^))+0{c-'). (4.8) 

This agrees with the definitions ( |3.1(]| ) and (|3.36|) of w|;(t, x) and W^(T, X) and the relation 



j:- = R-. (^a' -vi,a)+ 0{c~^), (4.9) 

which can be derived the definitions ( p.38|) , ( p.5|) -( pl6D of S", a^ and a with the help of the 
coordinate transformations ( p.44| )- ([3.45| ). 



2) External terms of order \'X.\^ require 

1 



B'it) = -^-vlv], + 2(1 + ^)w\^e) - (27 + 1) v],w{^e). (4.10) 



3) External terms of order |X|^ lead to 



1 



B''it) = -v%m. Q'^ + (1 + -i)w^''^\^E) - 7^g^''^ + ^7(5'^'#(x£;), (4.11) 



c^ R%k'j = -{l + ^)e,jkR\Ga 



-2(1 + 7) #*'^I(xe) + (1 + 27) i^i^ Wj](xe) + v'^R'^^^Qa + 0{c-'). (4.12) 



Eq. ( |4.12| ) relates the two functions Ga and i?°j. These functions together describe the spatial 



orientation of the local reference system with respect to the global one. According to ( [4. 12[ ) 
one can choose GaiT) so that i?"j = 5"^ and the resulting local reference system does not 
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rotate relative to the global one, i.e., the local reference system is kinematically non-rotating. 
Another possible choice C^, = results in a dynamically nonrotating local reference system, 
and the orthogonal matrix R"'^ in the transformations of the spatial coordinates represents 
the well-known de Sitter, Lense-Thirring and Thomas precessions. Note that ( [4.12| ) contains 



the well-known results for Lense-Thirring and de Sitter precessions in the PPN formalism 
(the second and third terms on the right-hand side, respectively). On the other hand, 
the fourth term represents Thomas precession appearing in special relativity and therefore, 
must be independent of the PPN parameters, which is indeed the case. The first terms in 
B^ and B^^ (Eqs. (|4.10| ) and (|4.11| )) come also from special relativity theory and are also 



independent of (3 and 7. 

4) External terms of order |X| , / > 2 give the external tidal vector potential in the local 
reference system 

W!^{T, X) = R% l^\t, x) - w\-K^) - w^i-x^^yi - v},W^{T, X) 

+ ^(^ {D'^'r^Erl - C,)] + 0{c-% (4.13) 



Eq. (|4.13|) together with (|3.15|) and (|3.43| ) leads to a Poisson-type equation for the unknown 
function C(t, x) 

C,u = 2b^''r% - 2(1 + 7) 44 = (7 - 2) 44- (4.14) 

This equation does not fix C(t, x) uniquely, but only up to a harmonic function. This fact 
reflects some ambiguity in constructing the local reference system in the form ( p.32|) -( ^.34|) . 



Assuming the function C together with its first and second partial derivatives with respect 
to X'^ to be equal to zero at X"- = 0, the general solution of (|4.14| ) reads 

C = Ci + C2, (4.15) 

where 

Ci(t,x) = ^4(d^ry (4.16) 
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is the particular solution of ([4.14|) , and 



^2(t,x) = E-[K....a,r^^..r^^ (4.17) 

p=3 "• 



Here V^^... a (t), p > 3 is a family of arbitrary STF tensors. Any closed-form solution of ( |4.14| ) 



leads to a closed-form metric tensor of the local reference system. E.g., we might choose the 
simplest solution C = Ci = (7 — 2) r| (a^r^)/10 corresponding to Ki...ap(^) = 0. In this 
case the tidal vector-potential reads 

W^{T, X) = R% \w\t, x) - w\^^) - wi.{^^)ri - v\, W^[T, X) 

27H^ i . fc _ ^7~^ , aVr|)+g'(c~'). (4.18) 

5(1 + 7) ^^ ^ ^^ 10(1 + 7) ^ ^J ^ ' ^ ' 



Note once again that the factor 7 + 1 in denominators of ( |4.13| ) and ( [4.18| ) cancels with the 



same one in the numerator of ( |3.32| ) and no rational functions of 7 appear in the metric 
tensor. 

It is easy to see from ( [4.13D -( P:.14| ) that independent of C2 the tidal vector potential Wf 



satisfies the equation 



ra / na -i i i''~\( „—'i\ 



W^,^, = T— ^ R% ^'e + 0{c-'). (4.19) 

1 + 7 



D. Matching of order 0{c ^) in qqq 

1) Matching of both terms C(c~^) and C(c~^) in ^fQO gives 

WeIt, X) + ^ *(r, X) = u.e((, x) (1 + 1 (7 + 1) 4) 

-H<l±i>u4(t,x).4 

2(/3-l 



C2 



w^{t, x) {w{^e) + a^ri,) + 0(c-^), (4.20) 



generalizing (f4.2|) to post-Newtonian accuracy. On the other hand, one can relate the local 
internal potential VFe with the global ones w-^ and w^ directly by using their definitions 
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( p.lOp and (|3.35|) . To this end, from the definitions ( p.37|) , ( |3.5| )- (|3l^ ) of E, a* and a, and 
the coordinate transformations ( p.44| )-( p^.45| ) one derives 



S(T, X) = a{t, x) (l + i ((1 + 7)4 + (2/3 - 87 - 1) (w(x^) + 4 4) 



--2(l + 7)t;i,a^ + 0(c- 



Using this the following result generalizing Theorem 5 of p2[ is found 



We{T, X) = WE{t, x) (^1 + i (7 + 1) 4) 



■^^^^wUt,^yE 



Here, 



+ - (2/3 - 7 - 1) W^it, X) (W(X£;) + 

-4(4/?-7-3)xJa^ + 0(c^^). 



X^(t, x) = i G / (T(t, x') |x - x'l rfV + 0{c~ 
2 jy 



a^r^ 



and 



X^(T, X) = ^G J^ S(T, X') |X - X'l rf^X' + 0{c-') = x^(t, x) + 0{c- 
A comparison of (|4.2CI| ) with (f4.22|) shows that 



v^(r, X) = -r^ {w^ (w(xs) + a^ r^,) - x!^ 4) + 0(c-') 



where 



r^ = 4/3 - 7 - 3 



(4.21) 



(4.22) 



(4.23) 



(4.24) 



(4.25) 



(4.26) 



is the well-known Nordtvedt parameter indicating a violation of the Strong Equivalence 
Principle in the PPN formalism. The fact that \E' is not equal to zero unless rj = means that 
requirement B from Introduction is violated and the internal gravitational field in our local 
reference system does not have the same form as the corresponding solution for an isolated 
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body. This has a profound physical meaning (see, also Section for further discussion) and 
is directly related with the violation of the equivalence principle for those theories where 
?7 7^ 0. In case of a system consisting of two spherically symmetrical nonrotating masses the 
term \l/ is equal exactly to the 6th and 7th terms in goo in Eq. (52) in |6^, providing an 
additional check of our results. 

Note that the expression ( [4.25| ) for \1/ results from the Newtonian potential Ue if we 



consider the expression ( [1.9|) for the "effective" gravitational constant as seen in the local 
system. 

2) Equating the external terms of order |X| give us the equation 

m = -^4 + 2(7 + 1) v'j,w\^e) -{l + D vIw{^e) + (/? - ^) ^'K)- (4-27) 

3) Matching the terms of order |X| leads to the equations of motion of the origin of the 
local system with respect to the barycentric PPN reference system 

a% = w,,{^e) - R^jQa (S^' - ^ {vl 5'^ + (2 + 7) wJ(x^) 5'^ + \v'ev'e)) 

+ ^(^2(1 + 7)#'(xi,) + (74 - 2(7 + I3)w{^e)) w,,{^e) - (27 + 1) v^eH^e) 

-2{l + ^)v'j,w%^E)-v],viw,{^E))+0{c'^). (4.28) 



Note that the Newtonian terms W/i_{^e) and Qa in (|4.28|) should be computed at moments of 
time t and T related to each other by T = t — 1/c^ A{t) + (9(c~^), that is by the coordinate 
transformation (|3.44|) evaluated at x = x^;. For Qa = 0, ( [4. 281) coincides with the PPN 



equations of motion of a test particle (geodetic equation), provided that the gravitational 
potentials of the local subsystem u7e, "U^e ^^^ excluded from the metric. 

We consider Qa as a known function of T specified in the local reference system (T, X°') 



with post-Newtonian accuracy. When computing the Newtonian terms in ( [4.28| ) with post- 
Newtonian accuracy, Qa should be expressed in terms of global coordinates. It is easy to 
see that all relativistic terms in ( |4.28| ) proportional to Qa are due to the re-computation of 



Qa from the local reference system to the global one. Indeed, if a test particle is situated 
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at the origin of the local reference system (Xfgg^. = 0) and if its local coordinate velocity is 
equal to zero {dX^^^^/dT = 0) then making use of ( p.44| )-( p.45| ) we get 



+0(c-^), (4.29) 

where d'^xl^^^/dt'^ — a^^ is the global coordinate acceleration of the test particle relative to 
the origin of the local reference system. The terms on the right-hand side of (|4.29| ) exactly 
match those in ( [4.28|) that are proportional to Qa- Qa is, thus, the acceleration of the 
instantaneous locally inertial reference system (whose origin coincides with that of the local 
reference system at a given moment of time) expressed in the local reference system. This 
also can be seen from the equations of motion of test particles relative to the local reference 
system (see, Section |VIII| below). The equations of motion can be represented as 



XLt = Q\T) + /"(Xtest, Xtest) + 0(c-^), (4.30) 

where /"(Xtest, Xtest) is some force depending on position and/or velocity of the particle (i.e., 
the only term in the equations of motion independent on the actual trajectory of the particle 
is Q°'). Note also that Qa can be directly related to the 4-acceleration of the worldline of 
the origin of the local reference system (see, e.g., p2| ). 

4) Matching of the external terms of order |X| , / > 2 in closed form leads to an expression 
for the relativistic tidal potential of the form 

W^{T, X) = wit, x) - w{-ke) - Wj(x£;)4 

+ i (^-2(1 + ^)v), {w\t, x) - W'(x£;) - tZ;:,(x^)r^) + (1 + 7)4W^t 

+ (1 + 7)^"'(x^) r^r^E + \iH^e) 4 + (^ - /? - 7) {a^ErE? 
+ (1 - 2/? - 27) QaX- ^eTe - 7 v),r'E ^,j (^e) 4 + ^7 4 i^'^.a^Q, 
+ C,t{T, X) + 2(1 - (3) (w(xs) + a^ r'^) Wt) + 0(0^^), (4.31) 

W^{T, X) = R\ l^w\t, x) - w'i^E) - w]^{^E)ri - vWt{T, X) 
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where C^t is a partial derivative of C with respect to the local coordinate time T. E.g., for 

C(T, X) = Ci(T, X) = ^X2 a^j, R- X^ (4.33) 



one gets 



C^t{T, X) = ^^X^ a% Rt X\ (4.34) 



Let us note that in the Newtonian terms in the right-hand side ( [4. 311) the moment of 



time t depends on both T and X according to ( p.44| )- (|3.45| ). This should be accounted for 
when computing, e.g., partial derivatives of Wt with respect to T. 



From ( [4.31| ) it follows that 



+ (7-1) (#(xs) + i?",a*^X")-2(/?-l)a| + 0(c-2), (4.35) 

W^T,aa66 = 0{c-^). (4.36) 



Eqs. ( |4.35| ) and ( |4.19| ) can be considered as effective field equations for the external potentials 



in local coordinates. 

We have derived all unknown functions with post-Newtonian accuracy and thereby ob- 
tained the local PPN metric tensor as well as the transformations between the local and 
global PPN reference systems. 

E. Energy-momentum tensor in the local reference system 



It is interesting to check the form of the energy-momentum tensor for an ideal fluid (|3.19|) 



in local coordinates, T""^. The energy-momentum tensor T°^ in the global PPN reference 
system is related to T°^ by 

^"'(^'X) = ^-^^T-(t,x). (4.37) 
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For an ideal fluid T^'^ is defined by (|3.19| ) or (|3.22|) . Substituting the transformation rules 
( p.44| )-( |3^ ) and the expressions (|i^), (|]6|) and (^?7D for A, D^^ and D^^'^ respectively one 
gets 



r°° = pc' 



i + ^n^ +n + 2i^ 



T^'^ =pcV'' 



1 + \[V^ + U + 2W 



+ 0{c- 



c 



a \rb I mb„ 



T"" = pV^V + d^'p + Oic 



^-2\ 



(4.38) 



where V"" is the velocity of matter relative to the local reference system, p{T, X) = p{t, x) is 
the invariant density, n(T, X) = n(t, x) is the specific energy density, and p{T, X) = p(t, x) 
is the isotropic pressure, W(T, X) = w(t, x) — u7(x£;) —a''^r^^ + 0{c~'^). One can also rewrite 
(g) as 



^00 ^ ^*^2 



1 + 



^ ^V2 + n-(37-2)W^') 



c2 V2 






r°" = p*cy" 



1 /I 



c2 V2 



l + ^[-V' + U-{3^-2)W 



+ 0{c~ 



+ -pV^ + 0(c- 
c 



^ab ^ p* ya yb ^ ^ab^ ^ 0{c''), 



(4.39) 



where 



p*(T,X)=p(T,X)(-^) 



1/2 '^"-^ 

(is 



1 + ^ f i ^' + 37W^ ) + C?(c 



(4.40) 



so that 



P^^A^-v')^o(.-^) 



(4.41) 



Thus, the functional form of the energy momentum tensor is the same in the global and 
local PPN reference systems. 

V. OTHER FORMS OF THE LOCAL METRIC 

In the previous Section we have succeeded to construct a local reference system (T, X") 
which satisfies property A from Section |, but does not satisfy property B. It is interesting to 
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check if there exist other versions of a local reference system satisfying properties A and/or 
B. To this end let us consider another set of local coordinates (TjX*^) related to (T, X") by 
the following transformations 

f = T, 

X" = X" + ^ ^ (r% R\ D'^X^ + R% R\ R\ D'^^X^X'') + C(c-^), (5.1) 

where D"^^ and D'^^^ are defined by ( |4.6| )-( |^I7| ). One can show that in {T^X"") the metric 
tensor reads 

Goo = -1 + 4 ^{f,iq -\f3w\f,iq + o(c-5), 



Goa = - ' ^3 " W\T,iq --T] (wi^E)X'' + X'^R''^d'j,X''--R\d'E^'j+Oic-'), 
Gab = Sab (l + 4 (7 W^(T,X) - ri (wi^E) + R'^.o'eX^))) + Oic-^), (5.2) 



W{T, X) = We{T, X) + Qa{T) X" + Wt{T, X), (5.3) 

VV'if, X) = W^{f, X) + ^ Babe C,{f) X^ + W^{f, X), (5.4) 



where 



We = G f f:(f,X') - ^ - d^X' 
Jv ^ ' IX-X'I 



/v "'" '""' IX-X' 

1 ^ 92 



IVe" = G j^ t'^if^X!) ^^_ d'X' + 0(c-2), (5.6) 

t = \ ("r"" + 7 T"'^ + 1 1°° (37 - 2/3 - 1) ly") + 0{c-^) , (5.7) 

E'^ = ir°'^ + 0(c-2), (5.8) 

c 

where T'^'^ are the components of the energy- momentum tensor in coordinates (T, X), and 
Wt and W^ are some external tidal (~ 0{X?)) potentials. 

The internal potentials We and We have the same form as the potentials of the cor- 
responding isolated source provided that the tidal influence of external metric is neglected 
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and, therefore, in coordinates (T, X) property B from Section | is satisfied. However, the 
terms proportional to r] in Goa and Gab represent non-tidal external potentials, thus violating 
property A. 

Although any reference system can be obviously used, we believe that the original local 
coordinates (T, X) are more convenient since a violation of the Strong Equivalence Principle 
is presented solely by the \E' term in Goo- Moreover, the appearance of non-tidal external 
potentials in Goa and Gab results in a more complicated interpretation of observational data 
in (T,X) (e.g., involving a complicated redefinition of locally measured physical constants 
and parameters). In the following we use a local reference system with coordinates (T, X). 

VI. MULTIPOLE EXPANSIONS OF THE INTERNAL POTENTIALS 

Although the metric tensor of the local PPN reference system described above enables 
us to discuss any physical phenomena in this reference system in the first post-Newtonian 
approximation, for many reasons an expansion of the gravitational potentials in terms of 
multipole moments is desirable. The main usefulness of mass (and spin) multipole moments 
of a body results from the fact that for astronomical, almost 'spheroidal' bodies usually 
only a few determine the equations of motion to sufficient accuracy. Numerical values of 
multipole moments can be obtained directly from various astronomical observations without 
evaluating volume integrals like ( p.35|) appearing in the local metric tensor. 

The complete solution of the problem of multipole expansions of gravitational potentials 
and equations of motion has been found in [^,^,^ for the case of General Relativity. The 
authors use the so called Blanchet-Damour (BD) mass multipole moments [^8|,^. The BD 
moments are defined physically by reference to multipole moments which can be read off 
from gravitational waves emitted by an isolated system of matter. The principal attractive 
property of the BD moments for their practical applications in multipole expansions is 



contained in Theorem 6 of [^]. The Theorem states that the multipole expansions of 
general-relativistic post-Newtonian gravitational potentials in terms of the BD multipole 
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moments are almost of Newtonian form p2|. In [^,^,^ it has been shown that this 
almost Newtonian form of multipole expansions is advantageous for both translational and 
rotational equations of motion. 

It turns out that in the PPN formalism we can find adequate multipole moments by 
"mimicking" the almost Newtonian form of multipole expansion of the post-Newtonian 
gravitational potential. It is easy to see that with the following definition of generalized 
BD moments 



M, 



V 2(2/ + 3) c2 dT^ Jv 



2(7 + l)(2/ + l) 1 d 



(/ + l)(2/ + 3) c2 dTJv 

Sl= I e^^^'^'X^-^^^T.^d'^X, I > 1 
Jv 



S"X"^d^X, / > 0, 



(6.1) 
(6.2) 



one gets 



°° (-1 
1=0 *■ 



MLdL^ + ^MLdL\X\ 



^'A^A,^Oic-% 



(6.3) 



°o f-iV 
2(1 + 7)W^^ = -2(1 + 7)GE 



1=1 



n 



■M-aL-ldL-lT^ + -r-—r£abc<ScL-ldbL-lT^i 
J^ L -\- L LA. 



-2(l + 7)A,a + 0(c-^), 



(6.4) 



where 



A-rV ^-^^' 2/ + 1 1 



(6.5) 



Vl= / T.^X'^^d^X. 
Jv 

The terms containing A can be eliminated from ( |6.3| 

coordinate 



(6.6) 
by a transformation of time 

(6.7) 



This gauge is called skeletonized harmonic gauge in |^ (no A-terms appear in 

and the gravitational potentials We and W^ are "skeletonized" by the moments Ml and 
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Sl)- The generalized BD multipole moments ( |6.1| )-( |OD obviously do not possess all nice 
properties of the BD moments in General Relativity. However, the expansion ( |6.3| ) has the 
same almost Newtonian form as the corresponding expansion in General Relativity. This 
fosters the hope that the generalized BD multipole moments are suitable for the use in 
multipole expanded PPN equations of motion. 



VII. TIDAL EXPANSIONS OF THE EXTERNAL POTENTIALS 



Although we have derived the local external potentials W^ and W^ in closed form, in 
many cases it is preferable to have tidal expansions of them. These expansions are very 
convenient, e.g., for a discussion of equations of motions in the local reference system. Tidal 
expansions of the external potentials entering Goo and Gqj are quite straightforward to derive 
and read 

oo-i -|0o-| -|00-| \^2 

QaX^ + W^ = Y.-Q,X'~-,T.n{RL-V,)x' + -j:--—K,X', (7.1) 



1=1 



^ 1=3 '" 



2(1 + 7)W^T = 2(1 + 7)E 



;e^ (/ + !)! 



£abc, ChL-l X 



OO -I 

+ T.niRaL-V,L)X^ 

1=2 ''■ 

+ ^(7-1) {R^a^j.X'^X' -2R\a^j,X') +0{c-'), (7.2) 

where the symmetric tracefree tidal moments Ql, Cl, Kl and Rl are defined as 

Ql = STFi jdi W(xs) + ^ (^2(1 + 7) R^i (Il-iw^-ke) - 2 (1 + 7) v),dLw\^E) 

+ (1 + 7) vl di w{xe) -T^lR^l R\ v'e 4 dbL-1 w{^e) + A^ 

2 

- 7 / + 2 (/3 - 1)] w{^e) dL w{^e) + [/ - 2 (1 + 7)] R^l v% dL-i W{^e) 

- 7/ (/ - 1) + 2 (/? - 1) / + 2 (1 + 7)] R'l a%dL^^w{^E) 
+ 0(c-^), />2, (7.3) 
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1 aft 



(2/3 + 27-l)i?«i?^a^ai;, 



A^ = 0, / > 3, 

i?i = 2 (1 + 7) STFi { i?"^ [rfi_i IIj'(xs) - v'e dL-i w{^e)] } + Oic-^] 

K={^- Ipi^E) + 2{l-P)al-7]Qa R% a% + 0{c-^) 
Ka = Qa -vQab R\ a% + {l ' I) R% a\ + 0(c-2), 
Kl = Ql-^ QaL R% a% + 0(c-'), / > 2, 



Cl — STFi l^ajbc-R'^i C^bL-l 



(wx 



X£;j -V^W 



(Xi.))}, / 



>2, 



(7.4) 

/ > 3, (7.5) 

(7.6) 
(7.7) 
(7.8) 

(7.9) 



(note that the expression in braces is already trace free). Here di is the shorthand notation 

Ql 



dLHt,^)=da,...a,Ht,^)=R%...R 



Ml 



X'^(^'X) 



(7.10) 



and dL^{yiE) should be understood as the corresponding derivative evaluated at x = x^;. In 



order to compute Ql for some given moment of local time T with the aid of ( [7.3| ) one has to 
evaluate d^ W(x£;) at moment of global time t related with T asT = t — 1/c^ A{t) + (9(c~^), 
that is by the transformations ( |3.44| ) evaluated at x = x^;. 

Expansion ( |7.1| ) with (|7.3|)-([7^) can be derived by differentiating ( |4.31| ) and using ( [4.351 )- 
( [4.36| ) and the fact that for any 5* 

Ql 



OO 1 

S* = V - " Q 

fr[ i\ 9X"i . . . x«' 
1 1 a' 

+ 77 



X' 



x=o 



S 



bb 



X^ X' + C(5 



bbcc) 



(7.11) 



x=o 



/! 2(2/ + 3) (9X«i ...X''^ 
The moments Kl represent the multipole expansion of (|4.35| ): c~^ Kl = W^^aaL- The 
quantities Vl are arbitrary STF moments appearing in the multipole expansion (|4.15| )- 
( [4.1?| ) of the function C{T, X). One can get rid of the moments Rl in ( |7.1|) - (|7.2|) by setting 
Vl = Rl, that is, by choosing a special solution for C(T, X) 



^ — 2 °° 1 ^ 

c(r, X) = -^x'R-, o^eX'^ + I: j,Rl X' 



(7.12) 
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In the limit of General Relativity /3 = 7 = 1 only the two families Ql and Cl play a role 
in (|7.1|) - (|7.2| ). These are related to the tidal moments Gl and H^ of the DSX formalism 

as 



Gl = Ql 



/3=7=1 



(7.13) 



Hr=-ACT 



13=1=1 



(7.14) 



Because of the gauge dependence the family Rl does not appear, e.g., in the equations 
of motion of test particles and A^ extended bodies which we derive in the next Sections. 
This can be shown also by computing the tidal expansions of the following quantities which 
directly enter the post-Newtonian equations of motion and are independent of the gauge of 
the time coordinate 



Q'^ + WT,a + ^2(l+7)iy^^ 



°° 1 ^1 °° 1 / 

E n QaLX"^ + - 2 (1 + 7) E 77 T-T ''^^^ ^^^-1 ^'''^' 
1=0 ''■ ^ 1=2 ''■''' ^ 

1-1 1 ^ ^2fL 1^1 (47 + 3)/-2(l + 7) ^ ^,^_i 



f^J\ 2(2/ + 3) 



^ 1=2 '■• 



2/ + 1 



+ ^ X« ( (7 - 1) #(x^) + 2 (1 - /3) a|) - ^ (7 - 1) i?"^, a^X^ 

2 - 1 °° 1 1 

+ — (7-l)i?^S^X'^''--ryE77 7T^^^^«^Q^^^ 



aL 



5c2 



c2 ' ^ Z! 2/ + 3 



^E 



1 1 



2c2 ' ;^ Z! 2/ + 3 



i?''.aEQa6LX^X^ + 0(c 



(7.15) 



-2(l + 7)£,,,iy^^, = 2(l + 7)E7y (aLX^ + ^^£afecX^''-'gcL-ij 



+ (7-l)£afeci?^ak^' + 0(c 



(7.16) 



In the limit of General Relativity ( [7.15| ) and ( |7.16| ) coincide with Equation (6.23a) and 
(6.23b) of @, respectively. 
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VIII. TRANSLATIONAL EQUATIONS OF MOTION OF A TEST PARTICLE IN 

THE LOCAL REFERENCE SYSTEM 

A. The equations of motion in closed form 

It is well know that a test particle moves along a geodetic of spacetime as determined 
by the metric tensor. The geodetic equation in the metric ( p.32| ) reads 

X: = W^a + ^{2(1 + 7) W^T - 2(/5 + l)WW,a - (27 + 1) W^tX^ - 2(1 + 7) W,tX',X: 
+7W^,aX,^X,^ + 2(l + 7) {w;i,-W',)x',} + 0{c~'). (8.1) 

Here, Xg are local coordinates of the test particle (e.g., a satellite), and all potentials and 
their derivatives should be evaluated at (T, X) = (T, Xs(T)). Accounting for the split 
( p.33| )-( |3.34| ) of the local potentials entering into ( p.32| ) it can be seen that the equations of 



motion in the local reference system read 

x: = $^ + $^, + 1 ($^-^^^p + $;^^ + <^l) + 0{c^'), (8.2) 

where ^^ ^^^ the terms describing the Earth's gravitational field in the absence of an external 
world 

$E = WE,a + 4 f 2 (1 + 7) WIt - 2 (7 + /5) We WE,a 



C2 



-(27 + 1) We,T X^ -2(1 + 7) WE,b X', X^ + 7 WE,a X's ^s 

+ 2 (1 + 7) X', {Wl, - Wi:)) + 0{c~'), (8.3) 



The Earth-third body coupling term $"oup due to the only nonlinear term in (|8.1| ) reads 

•^coup = -2 (7 + /?) (We {Qa + WT,a) + W^E,a {QaX^ + W^t)) + 0{c-^). (8.4) 

The "gravito-electric" part <l>gi of the influence of external bodies (independent of the velocity 
of the satellite) reads 

$:i = Qa + WT,a + 4 (2(1 + 7) WIt - 2(/3 + 7) {Qb X'^ + Wt) {Qa + W^T,a) 



C 

-4 



+il+^)eabcax:]+Oic-~^), (8.5) 
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and the "gravito-magnetic" part $|^g (depending on the velocity of the sateUite) of the 
influence of external bodies is given by 

$:^g = - (1 + 27) {wt,t + QbX^) x: -2(1+7) (w^T,6 + Qb) x'^x: 

+ 1 (W^T,a + Qa) Xl Xl + 4(1 + 7) lyj?''] Xl + 2(1 + 7) Sabc C, X^ + 0{c-^), (8.6) 

Let us remind that $"1 and $^ contain also the two arbitrary functions Qa and Ca that 
characterize the world line of the local reference system and rotation of its spatial axes. 
A reasonable choice of Qa leading to vanishing of the BD mass dipole A^" of the central 



body will be discussed in Section IX E . As it was noted above (see the discussion after 



( 4.12| )) a dynamically nonrotating local reference system corresponds to Ca = while in a 



kinematically nonrotating local reference system Ca contains the geodetic, Lense-Thirring 
and Thomas precessions. 

The term $|, represents the additional acceleration due to the term (^4.25 ) 



$^(T,x) = ^,a(r,x) 

= -V {W^,a {w{^e) + a^E Te) - X%{T, X) R\ d^ + W^ R% a*^) + 0{c-^). (8.7) 

In case of two spherically symmetrical nonrotating masses the term \1' (the potential of $|,) 
coincides with the 6th and 7th terms in g^Q of Eq. (52) in |6^. A discussion of observability 



of this effect in satellite motion data can be found in Section VI of 34 . 



B. Multipole expansions of the equations of motion 

For many applications it is convenient to expand the equations of motion derived in the 
previous Section into multipole moment of the gravitational field of the central body and into 
tidal moments of the external potentials. E.g., it is a common practice to use the equations 
of motion of near-Earth satellites in expanded form and to fit the multipole moments of the 
Earth's gravitational field. Substituting the expansions derived in Sections |VI| and |VII| into 
and ( |8.7| ) one gets 
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°° f-lV 1 



n 
1=0 '■• 

1 oo 



X.. 






2i^ + P)MLdaL 



1 



Xs 



V fc=0 



Xsl 



— . „ — ^abc <SbL dcL 



k\ 



MkOk 



1 



Xs 



+ ^?^ X'^ S,L Bdcla db],L ^] + 0{C~^' 



1 + 2 



|X« 



XI 



(8.8) 



$" 



lY 1 



coup 



-2(7 + /3)G EE^7^i:T-^LQai.Xf9z. 



/! A;! 



vi=0 fc=0 

oo oo /' -| \l -| 



x« 



+EE 

Z=0 fc=l 



/! A;! 



-A<iQ;,Xf 5,L 



IX. 



+ 0{c~'), 



(8.9) 



1 / 



■-s 



1=0 '• '^ \ z=i '• ' + -■- 

^1 1 /s ^2tL ^1 (47 + 3)/-2(l + 7) ^ ^ 

6^ /! 2(2rr3) ^-^^ ^^ ^^ "^/! 2m ^^-^ • 



aL-l 

s 



Z=0 

+^X:^i^-l)tji^E) + 2{l-P)al 
2 - °° 1 1 

1 °° 1 1 

"2''§/!27T3^''^"^^"''^^'^^ 

oo oo 1 -| 

-2(7 + /5)EEn7:yQ«^Qi^^s^^. 
i=0 /c=l ''■ '^■ 



6 



(7-l)i?",a'sX 



-R j a^ QbL ^; 



s 



2 i>L 

s 



s 



0{c- 



^.10) 



$ 



mg 



oo 1 oo 1 oo 1 



1=1 



l\ 



1=0 



l\ 



1=0 



°° 1 / - 1 -, . 

-2(1 + 7)5.,,X,''^- (c^LX^ + j^e^aeXt'^Q^L 



a • j \rb vb\ 



+ (7-1) {X^ R\ d^E ^s - R% a'E ^s ^s ) + 0{c 



^-2\ 



(8.11) 



where we use ( [7.15|) -( [7rT6D and the identity 2 W+'' = Sabc^cde W- 



a,b] 



id,e 
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The multipole expansion of $|, reads 

CO / 1 y f / 1 



1=0 '■• *- ^ l^sl 

+ -R"i (^E ^L prrr — 7:R i<^^E ^abL |Xs 

1 .,„.,„ 1 



2(2/ + 3) 
where 



AfLR\a%dabL T^ri] + C?(c-'), (8.12) 



7Vl= / SX^X^rf^X (8.13) 



V 



Here we use the multipole expansion of x {T, X 



1 ^^(-1)' / . . . ,,,, 1 ., . 1 



X^(T,X) = -G5:^ (■^^^^l^l + 27T3-^^^^pqj+^^'~'^' ^^-^^^ 



which can be verified directly from ( |4.24|) . The following relations are useful for evaluation 
of the expansions derived above 

^^^ = (-!)' (2^-1)"^' (8-15) 

d, |X| = i-iy-' (2/ - 3)!! ^ (^N^ - ^^^ [^[^^^^_ ^^ <5("-^''^iV^-2)^ , (8.16) 

where N"" = X''/\X\ and N^ = N^^'^N^^ . . . X"'>. The shorthand notations (9l(1/|Xs|) and 
(92,|Xs| mean (9^,(1/ |X|) and (92,|X|, correspondingly, evaluated at X = Xg. 

It is important to note that if r/ 7^ the equations of motion of a test particle contain 
not only the BD moments A4l and Sl as in case of General Relativity, but also one more 
family of mass moments A/l defined by ( ^.13|) . 



IX. TRANSLATIONAL AND ROTATIONAL EQUATIONS OF MOTION OF N 

EXTENDED BODIES 

Translational equations of motion of the origin of the local reference system is given by 
( [4.28| ). In order to relate the origin of the local reference system to the post-Newtonian 
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center of mass of the body (i.e., to fix Qa) we can employ the local equations of motion in 
the local PPN reference system 



-a/3 



r°% = 0. (9.1) 



A. Local equations of motion: T ^ = 

Let us consider the integral 

J {-g)T'^^d'X = 0. (9.2) 

In order to evaluate ( |9.2|) we have to calculate the Christoffel symbols as well as the metric 
determinant for the local PPN metric 

n, = --^W, + 0{c-'), (9.3) 

Tl = -^^W,^ + 0{c-% (9.4) 

C = ^ ((1 + 7) {w], + W':}+j6-'W,t)+0{c-'), (9.5) 

-g = l + \2{3^-l)W + 0{c-^). (9.6) 



Substituting (|]3D-(|]6D into (|]|) one gets 

-^Mt = F + 0{c-^), (9.7) 

where M^ is the Tolman or ADM mass 

Mt{T)= I p^jT,X)d'X, (9.8) 

with 

1 //I 



P,JT,X) = h + _(^(^- + 3(7 - 1)) We 



+ (3 7-2) (g, X^ + Wt)]]^ r°o + 0{c~'). (9.9) 



In case of an ideal fluid 
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p^jT,X) = p* [1 + ^I^^V' + U-^We)) +0{c-'), (9.10) 



p* being defined by ( [4.40| ). This justifies that p^^j as given by ( p.9| ) is indeed the usual 



Tolman density. However, the definition ( p.9| ) is independent of the model for T^'^. The 
right-hand side of (|9.7|) reads 

F=\ I T.'' {Qa + W-T^a) d^X + 0{c-^). (9.11) 

c^ Jv 



For an isolated body the mass Mt is obviously conserved. The BD mass Ai defined by ( |6.1|) 
is related with the Tolman mass Mt by 

Al^M. + l,(^. + li(,-l)|,|sxVX 

+ 4(1-7-2/3) Qa / SX^rf^x 
c^ Jv 

+ \ / S ((1 - 2/5) VTt - 7 WT,a X'') £X + 0(c-^), (9.12) 

c^ Jv 



where 



n^ = -]- I T.WEd^'X + Oic-^) (9.13) 

2 Jv 

is the internal gravitational energy of the body. One might have expected the appearance 
of the second term (containing ^2^;) in (|9.12|) since Mt is an inertial mass and M. is the 



active gravitational mass appearing in the expansion of the gravitational potential ( |6.3|) of 
the body. The third term can be related with a Newtonian virial theorem. It vanishes for 
an isolated body under the assumption of secular stationarity. The last two terms in (|9.12|) 



are proportional to Qa and the external tidal potential P^t respectively, and vanish for an 
isolated body. We see that the local equations of motion do not constrain the BD mass 
M. as was the case in General Relativity, but only the Tolman mass Mt- Moreover, in the 
framework of the PPN formalism the BD mass M. is not conserved even for an isolated body 
unless some specific physical conditions are met (secular stationarity and f^E = 0). 
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B. Local equations of motion: To = 



Let us now consider the integral 



V 



-g)f''%£x = o. 



(9.14) 



The Christoffel symbols necessary for the evaluation of ( |9.14| ) read 



Ko = -\ Wa + ^(iP + ^)WW,a-il + 7) W^^ + Oic-'), 



9 i" ' A 



Kt = 4 ((1 + ^) {K - ^i) + 7 5^' Wt) + 0{c 



ab 



r^ = ^ 7 -5 w,a + (J"" W,, + 5^' WA + 0{c 



(9.15) 
(9.16) 
(9.17) 



Substituting (p[5D-(^T^) and (|J) into [ W^M one gets 






M^ = F^ + -F^ + 0{c-''), 



(9.18) 



where 



J y 



(9.19) 



Jv 



d'X + \^ f g'^d'X + 0{c-'), 
c^ dT Jv 



(9.20) 



r = j: {Q^ + WT,a) + ^ 2(1 + 7) (^S Q 6abc Cb X' + W^^ 



+£"(-£,,, c, + iy^,,-iy^J), 



^» = S" (Q, + l^T.b) X'' - (27 + 1) S'^ {Q, X' + W^t), 



(9.21) 
(9.22) 






(9.23) 



It is easy to see from ( ^.18D -( p.23| ) that for isolated body M^ vanishes. The BD mass 
dipole A^" defined by (|0| )and M^ are related by 
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where 



+^ (1 - 2/?) / s (g.x'^ + vTt) x"^ (i^^ + 

2 Jy 



(9.24) 



(9.25) 



If ?7 7^ or 7 7^ 1, the second time derivative of the BD mass dipole M!^ does not vanish 
even for one isolated body unless specific physical conditions are met (secular stationarity 
and ii|). 

C. Local equations of motion: Eahc^^T n = 

The local PPN reference system of an extended massive body described above allows 
us to derive also rotational equations of motion of the body in the framework of the PPN 
formalism. Substituting ( p.l5| )-( |9.17| ) and ( p.6| ) into 



eabcJ{-G)X'r''%d'X = 



(9.26) 



one gets the following rotational equations of motion (see, p8| , |70| for a more detailed dis- 
cussion) 



A.S^ = L- + ^^L% + 0{c-') 



(9.27) 



where the PPN spin S"' is defined by 



S'' = eabc f X'p'd^X + Oic- 



(9.28) 



with 
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p'^ = SMl + ^^iy)--^GS / S^(T,X') ^^^ ^^^ '^Vt ''"'' cf'^' + 0(c-^), (9.29) 
c^ 2c^ jy |X — X I 

Here, 

L» = Sake I X' f (fX + 0{c-^), (9.31) 

and 

L%=eab, f i:x'^,,£X = T]eabAR"ia'E + Oic-^), (9.32) 

J Y 



where f^ and Q% are defined by ( |9.21| ) and (|9.25| ), respectively. For one isolated body S"" is 



conserved. It is easy to see that S"" coincide with the BD spin moment S"" defined by ( p.2| ) 
in the Newtonian approximation 

S'' = S'' + 0{c~^). (9.33) 

The torque L^ comes from the term \E' in the local PPN metric and represents an analogy 
of the Nordtvedt effect in the rotational equations of motion of an extended body in the local 
PPN metric. These rotational equations of motion and the torque L% has been discussed in 



a previous paper p^ . 



D. Multipole expansions of the relations between the Tolman and BD masses and 

mass dipoles 

For further analysis it is interesting to expand explicitly the relations for the Tolman 
mass Mj- and mass dipole M^ with the BD mass Ai and mass dipole A^" defined by ( |9.12[ ) 
and ( 9.24 ) respectively. One gets 



1 °° 1 
+ - E 77 (1 - 2/3 - 7 -Ml Ql + 0{c-'), (9.34) 



c2 "^ /' 
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where 



1 °° 1 

+ -E n (l-2/?-70-^aLQL + 
C ^^-^ i. 



Va= / S''X"^rf=^X (9.36) 

JV 



Here, A/l is defined by (|8.13| ). The mass and mass dipole are the only moments which 



appear in two different forms in the formahsm: 1) the BD mass and BD mass dipole in 
the multipole expansions of the internal gravitational potential (|6.3| ), and 2) the Tolman 
mass and Tolman mass dipole in the left-hand side of the local equations of motions ( p.7| ) 
and ( |9.18| ). The physical significance of the BD moments lies in their appearance in the 
multipole expansion of the post-Newtonian gravitational field, that of the Tolman mass and 
mass dipole lies in representing integrals of local equations of motion for one isolated body. 
In General Relativity the BD and Tolman mass and mass dipoles coincide for isolated bodies. 
In contrast to this, within the PPN formalism this is no longer the case and the BD mass M. 
and mass dipole A^" and the Tolman mass M^ and mass dipole M^ are essentially different 
objects. Higher multipole moments appear only in the multipole expansions of the internal 
gravitational potential (|6.3|) . 

E. Multipole-expanded local translational equations of motion 

Using the definitions of the BD multipole moments and external tidal moments intro- 



duced in Sections ^ and [VI 1| one can derive the following multipole expansions 



M = AM 



1 OO 1 



C i^^ i. 
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+ (/ + (7 - 1) / + 2 (/3 - 1) MlQl + 0{c~^), 



(9.37) 



Ma = AMa + R% a^E ^M 
1 

n 



oo 1 

J2n^LQaL 
1=0 



1 °° 1 / 1 1 

— 2 (1 + 7) ^ - ( -^^ Sl+1 CaL+l + ^-^ Babe MbL CcL 
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n V/ + 2 
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£abc-MbL CcL 
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£abc SbL QcL 
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(/ + l)(2/ + 3) 



£abc SbL QcL 
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- E 77 I ^^^ + (7 - 1) ^^^ + 2 (/5 - 1) ) A^aL Ql 



c- 



/! 



1=1 

1^1 1 



/ + 1 



/ + 1 



"c^"'^/! 2/ + 3 



MaLR^a^EQbL 



1 °° 1 1 
-2?^g/!27T^^'^«-^^^"^^ 

2c2 "^ ^Q /! 2/ + 3 rfT2 1,^^^"^ 
1 , rf2 

6c2 ^ ^ ^ rfT2 

1 



ATi^'^.a^ 



3c2 



+ T^Ma (7-1) f/(t,xs) - 2 (/? - 1) a^ + — (7 - 1) i?", a^^AI 



d KAOb 



5c2 



+ - (7 - 1) - A^,5 ^ i Ci^E - o ^<^bc '5^ R"i «^ + (^(C 



(9.38) 



with 



oo 1 

'^'^ = E 77 ^«fec MbL QcL + C'(c 



-2> 



(9.39) 



and 



AM 



r^^E + -{l-l)M 



(9.40) 



AM, = 1 fr^f]^ - 1 (7 - 1) AT, + I (7 - 1)P„ 



10 



5 



(9.41) 
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AM and AM^ represent the internal parts of the difference between the gravitational BD 
mass and mass dipole, A4 and A^", and the Tolman mass M^ and M^ given by ( |9.34| )- 
(g). In the limit of General Relativity (|9:37|)-(|09D coincide with Eqs. (4.20a)-(4.21c) 
of[0. 

Analogous to the translational equations of motion discussed in Section |VTT|, Eqs. ( |9.37| )- 



( |9.38| ) contain not only the two families of BD moments Aii and Sl, but one more family of 
multipole moments A/l. Moreover, (|9.37|) - (|9.38|) contain also a number of additional terms 



that vanish in General Relativity, but do not vanish for rj = 0. The physical meaning of 



some of these terms will be discussed in Section IX G below. 



F. Multipole-expanded rotational equations of motion 



In the same way for the post-Newtonian rotational equations of motion one gets 

°° 1 / 1 / + 1 \ 

S" = Babe E ^ [-^bL QcL - ^ 2(1 + 7) ^ SbL CA 

with 

^A = ^* - 4^ (7 - 1) eabc K R% a], + 1(^-1) Sabc n R% a% = c^ Eabc AMfe R\ a%, (9.43) 
lU 5 

SI = -2(1 +^)Y^--^—MaL Cl, (9.44) 



^-2^^ + ^)g^27T^^^^- 

°° 1 (21 + 3) 
+2(1 + 7) Sabc E ^ (; + 2)(2/ + 5) ^^^ ^^^ 

^1 / + 2(7 + 2) / . ■ ■. ^ 

^"'^^^ g n 2(Z + 2)(2/ + 5) V-^^" ^^^ " -^'^ ^^^ 

°° 1 1 (i 

+ (1 + 7) Sabc E 77 // , o^Mr , ^^ 1^ (^^l Qcl) 



^Z! (/ + 2)(2/ + 5) dT 
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!<- 


^)£abc'PbR'^iO.E 






^- 


l)eabc{MkR\a%- 


-KR\ 


a^E 




-l)e.,e^ {U,R\ 


o^e) 





(9.45) 



and 



°° 1 1 
L^, = -V Bate E yy 2(2/ + 5) -^^^ ^'^ ^^ ^^'^^^ ^9-46) 

Here, L^ is defined by (|9.32|) and Vl by ( |6.6|) . We prefer liere to consider L^ as part of 
L'^ which is proportional to the difference of the BD and Tolman dipoles AM^. This might 
provide a deeper insight into the nature of this torque. As in General Relativity we can now 
change the definition of the post-Newtonian spin into 5"" = S"" + S^ + S2 to get rid of S^ 
and 5*2 on the right-hand side of ( |9.42D (see, |^ for a discussion of this point): 

00 1/ 1 ] 4- ^ \i 

S''' = SabcT^J^i-^bLQcL- ^2{l + j) j^S.lC^l) + ^ {1% + L^) + 0{c"') . (9.47) 

The torque L° is in some sense similar to L%: it is proportional to the Nordtvedt parameter r) 
and the acceleration of the body relative to the global reference system a^. For that reason, 
it is related with a violation of the Strong Equivalence Principle. L% has been discussed in 
detail in a previous paper [^ . 



Let us note here that the rotational equations of motion and the definition of the post- 
Newtonian spin should be considered as being formal in the first place. Further analyses 
will be needed to see how these concepts can be used efficiently in practice. 

G. Mass monopole-spin dipole approximation 

Let us discuss the skeletonized equations of motion of A^ extended massive bodies with 
full multipole structure derived above for the case that the A^ bodies possess only mass and 
spin. The corresponding mass monopole-spin dipole model for the multipole structure of 
each body B can be mathematically formulated as 
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Ml = 0(c-^), 

Sl = 0{c-'), 

AfL = Oic-'), 



I > 1, 
/ > 2, 
/ > 1, 



(9.48) 
(9.49) 
(9.50) 
(9.51) 
(9.52) 



These assumptions agree with what we expect from a naive Newtonian model of spherically 
symmetric rotating body. From (pD-^H), (lOTD-dOsD, (|^-(|^) and (^AEjj-j^I^ 
one finds 



AMR^y^ + MQa --^vAfR\o'EQab + ^ (1 - 7) ^ (A/'^"i4 
-\ (1 + j)St,Cab -Mil- l)eabcSbR%a^E = 0{c-'). 



(9.53) 



(9.54) 



(9.55) 



(9.56) 



c^ ■ • c^ 2 

Eq. ( p.56| ) allows us to derive Qa for the case that the origin of the local reference system 
coincides with the BD center of mass of the central body [Aia = 0). For convenience we 
split Qa into a part independent of the spin of the body Q™ (that is, Q™ corresponds to 
mass monopole model for the central body) and a part proportional to the spin of the body 

Ql 

Qa = Q7 + Ql (9.57) 



^« ~ 'M' Ri^E + -^Vj^RiaEQab 



^(l-7)i:l,(A/'i?" 



6c^ 



M dT^ 



i"'E 



0[c 



(9.58) 



Ql = \^{X^l)^S,C,,-^^{\-^)^eabcS,R\d\^0{c'') 



(9.59) 
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From assumptions (|9.48|) - (|9.49|) and expansions ( |6.3|) -( |674D in the local PPN reference system 
of each body B we have 

Wb = GMbt^ + 0{c-'), (9.60) 

\^b\ 

W^ = lGeabcSI,dbT^ + 0{c-'). (9.61) 



Here, for each body the skeletonized-harmonic gauge (|6.7| ) was used (see, [^ for more 
detail). X^ are spatial coordinates of the local PPN reference system of body B. Now we 
derive the potentials of body B in the global PPN metric following the steps of Appendix 
C of m. From j ^AEj ), (pj) and g^) one gets 



- fl - i fi (x],n:^' + ^Wsi^B) + Iix^bTb ) ) + (^(c-'), (9.62) 



|Xb| tb V 6 

where r^ = x* — a;^, n^ = v^b/^b- Here, x^ are the coordinates of body B in the global PPN 
reference system, and wb is the external potential appearing, e.g., in (|3.9|) when constructing 
the local PPN reference system for body B. Then substituting (|9.60|) - (|9.61|) and ( |9.62| ) into 



( [4.20| ) and ( [4.8|) we get the explicit form of the gravitational potential of the body B in the 
global PPN reference system 

WB = ^ (l + ^ ( (^ + ^) *^ - ^2/3 - 1) Wb{^b) 

+ ^ [^B ~ ^b'^'b ~ [^B^b) 

+4 (1 + l)Ge,,kx'B 4 dj — + 0(c~^). (9.63) 

C IB 

Here, s^b — ^i ^b ^^^ ^i is the rotational matrix i?" appearing in the local PPN 
reference system for body B. The expression for the vector potential W^ in the global 
reference system can be derived immediately from (|4.8|) , ( |9.63| ) and ( p.61|) 



w'b = ^^x^B + lGei,ks'hd, — + 0(c-2). (9.64) 

Tb 2 rB 

Now we note that in analogy to ( ^.9] ) the gravitational potentials of the global PPN metric 
can be split into the sum of the contributions of the bodies 
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B 


B 



(9.65) 

Hence we get the explicit form of the gravitational potentials of the global PPN reference 
system within our mass monopole-spin dipole model ( |9.60| )-( p.61| ) 

^GMb (, iff 1\ .2 /oo ^, ^ GMc 

E-— 1 + 3 7 + 77 4 -(2/3-1) E 



B ^B \ c y\ 2/ ^_^g rcB 

+-^l^T.GMBrB,tt 

+ 4 (1 + 7) T.Ge^jkx'Bs'hd, — + 0{c-'), (9.66) 

w' = E^^4 + I T.Gsi.ks'hd, - + 0{c-'), (9.67) 

B ^B 2 B ^B 

where tcb = xc — x^ and rcB = ItcbI for any C and B. 

The equations of translational motion of body E ([4.28|) can be now written as 

4 = «km + 4,s- (9-68) 

Here the part a^ j^^ depends only upon the masses of the bodies 

(^E,ra = - E ^ -^B — -R" Qm 

Bj^E '''eB 

2 



3(4^ 1^^^ 






^ ' EB ^ C^E.B ' BC 

- (1 + 7) x-^ x-^ - 7 i:-^ x'b + 2(1 + 7) x;^ i;-^ ^ 
+4 EG-^B%|2(l + 7)ii-(27 + l)ii.|(il.-iy 

^ B^E ^EB ^ J 

4 (2^ + I) ^ T^ E G A<o ^ + O(c-). (9.69) 

These are the PPN equations of motion for a system of mass monopoles. They agree with 
the well-known PPN Lorentz-Droste-Einstein-Infeld- Hoffmann equations of motion (see, e.g., 
|1[]). Finally, a^g are additional acceleration terms due to the spins of the bodies 
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(^E.s - --R^i Qs 



(1 + 7) G Y. Uijk 4 (i^ - is) dji + Sjik 4 (i^ - ifi) ^i 



C ^^^ \ Teb fEB 

+0{c-^). (9.70) 

The geodetic deviation Q1 due to spin of body E can be further simphfied with the mass 
monopole-spin dipole model we consider. From (|7^) , (|9.66| ) and ( |9.67|) one gets 



^ iQs ~ 2 A ^ ^^i'^ ^E ^E 
C JVIe 



C JVIe b^E ^ ^EB 

\ £ijk Se \^E -^b) ^jl 

fEB 

+^ (1 + 7) ^ G 5: s^E s'b d^,u -^ + 0{c-'). (9.71) 

IC JVIe q_^^ Teb 

Here we used the identity eij[k Aq^ = | ^kim Aij valid of any trace-free Aij (so that, An = 0). 
Note that 

a.^^ '^^^^^y''"^" , (9.72) 

I EB I EB 

Q ^ 1 _ 3 r|^ {t'es S^'' + r^EB ^"^ + ^Ij? ^'^) - 15 r%B rU ^eb .g -^gN 

'^'^ rEB rls 

Equations (|9H)-(|97|) together with j ^M ) for 7 = /3 = 1 agree with Eqs. (6.29)-(6.36) 



of [^. These equations of motion can be derived from Lagrangian quoted in Section |. 



Note that we did not assume the constancy of the Blanchet-Damour masses Mb- Strictly 
speaking one has to account for the time dependence of the masses A4b during a [numerical] 
integration of the equations of motion ( |9.69| ). It is also clear that this effect is very small 
and usually may be neglected in practical applications. Another point is that the Nordtvedt 
effect here is represented here not only by the terms rjQB and rj^E, but through the whole 
difference of the BD active gravitational mass Ai and the inertial Tolman mass Mt of each 
body. According to ( |9.54D and ( |9.4UD the difference contains not only t]Q, but also terms 



(7 — 1)N'/Q that do not vanish even for spherical bodies. 
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Taking into account that the BD mass Ais for each body is the active gravitational mass 
(see, ( p.l2| )), it is easy to show that these equations of motion coincide with Eqs. (6.31)- 
(6.34) and (6.47) of |I| if only the two PPN parameters (3 and 7 are retained there, all 
total derivatives of any internal integrals are dropped in our equations (secular stationarity 
is assumed in []I|) and the A/'-dependent terms are neglected in our equations. The M- 
dependent terms proportional to t] and to 7 — 1 can be explained by a coupling of the body's 
extension to the background Ricci tensor, as already mentioned in Section | and coincides 
with the corresponding terms derived by Nordtvedt in ITIHT^,^ . 



The spin-dependent part of the acceleration ( p.70| )- (|9.71| ) coincides with the known spin 



orbit coupling derived, e.g., in []60| , |6l| , |53| and can be derived from Lagrangian ( |1 . 1 7|) . The 



first term in ( p.71D is independent of the PPN parameters. It results from the well-known 



Thomas precession and has been first discussed by Damour [Q (see, also [|6T|). 
For the mass monopole-spin dipole model (|9.48| )- (|9.52| ) extended by assuming 



Vl = C(c-'), / > 2 (9.74) 

the rotational equations of motion ( |9.42| ) read 

S"" = -^ (1 + 7) Babe '^^ ^^ + ^ ^ (1 + ^) ^ (^^'^) + ^(^"')- (9-75) 

The second term in the right-hand side of (|9.75|) can be absorbed by changing the definition 
of the post-Newtonian spin as it was done in ||2^ for General Relativity. 

If the local reference system is chosen to be kinematically nonrotating, R^ = 6^, then the 
first term in (|9.75| ) gives precession of the spin due to geodetic, Lense-Thirring and Thomas 
precessions (recall that in a kinematically nonrotating local reference system R'^^ = 5°"^ and 
Ca represent just the rotation of the local reference system relative to the global one due 
to geodetic precession (see, ( [4.12| )). If, however, a local dynamically nonrotating reference 



system is chosen then C^ = and the term vanishes. Thus, we prove that in the PPN 
formalism an extended massive body with mass and spin only undergoes the same relativistic 
precession in the background gravitational field as test massless gyroscope. This generalizes 
the corresponding result for General Relativity given in |2^ 
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X. OUTLOOK AND CONCLUSIONS 

In this paper we have formulated a new framework of relativistic celestial mechanics in 
the first post-Newtonian approximation with PPN parameters 7 and (3. Here we list some 
of the most interesting results that were derived in the main part of the paper. 

• It is impossible to construct a local reference system satisfying simultaneously proper- 
ties A and B formulated in Section | unless r/ = 0. Either property A or property B 
separately can be satisfied. For practical applications the local reference system where 
the external gravitational potential is represented by tidal effects only (i.e., property 
A is satisfied) seems to be more useful. 

• As in General Relativity it is still possible to formulate the theory in terms of two 
material variables a and a* related with T"^ by 



It is possible to find an empirical definition of mass multipole moments (|6.1| )- ([6l^ ) that 
allows us to keep the same almost Newtonian multipole expansion ( |6.3| ) of the local 
gravitational potential W as in General Relativity. Those multipole moments gener- 
alize the Blanchet-Damour moments introduced in [^ in the framework of General 
Relativity. 

Unlike in General Relativity the BD-like mass and dipole moments (gravitational mass 
and mass dipole) do not coincide with the Tolman (or ADM) mass and dipole moment 
(inertial mass and mass dipole) even for one isolated body, where the BD-like mass in 
general is not conserved. 

The tidal expansions of external potentials in the local reference system involve (3 
and 7. All tidal terms can be divided into two groups: those which exist in General 
Relativity (coefficients now are polynomials of (3 and 7 rather than rational numbers) 
and those vanishing in General Relativity. 
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Equations of motion of a test particle in the local reference system contain terms 
proportional to the acceleration of the central body relative to the global PPN reference 
system. The effect is proportional to rj and is a consequence of a violation of the Strong 
Equivalence Principle. 

For r] ^ the equations of motion of a test particle in the local reference system contain 
an additional family of multipole moments A/i, defined by (|8.13|) . Those moments 
also characterize the gravitational field of the central body, but do not appear in the 
multipole expansions ( |6.3| )-( |OD of post-Newtonian gravitational potentials W and 

Wa. 

For P and 7 unequal to 1 equations of motion of both test particles (satellites) and 
massive bodies involve many additional terms proportional to r], 7 — 1 and /? — 1 which 
vanish in General Relativity. The terms proportional to rj come from a violation of the 
Strong Equivalence Principle, some other terms result from a coupling to the external 
Ricci tensor. 

It is no longer possible to construct a "point-mass" approximation to the equations 
of motion of A^ extended bodies by assuming that all the bodies have constant BD 
mass and all other BD mass and spin multipole moments vanish identically as was the 
case in General Relativity [^. A number of additional assumptions concerning the 
moments A/l as well as Va and Qa is indispensable. 

If f3 and 7 are unequal to 1 local equations of motions cannot be represented as a 
bilinear series of BD-like mass and spin multipole and tidal moments: one more family 
of multipole moments Afi and a number of additional terms (containing Va, ^ and 
fi") appear here. 

It is possible to derive equations of motion of A^ massive extended bodies possessing 
only mass monopoles Ai (and moments of inertia A/"), and, possibly, spin dipoles iS\ 
without assuming secular stationarity. These equations of motion do not assume the 
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gravitational masses of the bodies to be constant. The equations coincide with the 
known results (e.g., the EIH equations of motion) in the corresponding limits. 

Massive bodies possessing only mass monopoles Ai (and moments of inertia Af) 
no longer move along a geodesic of the background metric. This is usually called 
Nordtvedt effect and presents a direct consequence of a violation of the Strong Equiv- 
alence Principle. As compared with standard treatment of the Nordtvedt effect [|l|, 
the geodetic deviation of such a body contains not only terms rjQ proportional to 
internal gravitational energy of the bodies, but the complete differences between the 
gravitational and inertial masses of the bodies (|9.4(J|). 



Rotational equations of motion contain a term analogous to the Nordtvedt term in the 
translational equations of motion (i.e., a term proportional to t] and depending on the 
acceleration of the body relative to the global PPN metric). This question is discussed 



in 1661 in more detail. 



The framework presented here extends the Brumberg-Kopeikin and DSX- formalisms by 
the introduction of PPN-parameters. It also extends the classical Nordtvedt- Will PPN- 
framework by: the formulation of a theory of astronomical reference frames with new trans- 
formation rules between global and local coordinates comoving with a body of the grav- 
itational A^-body system, improved definitions of mass and spin multipole moments that 
characterize the gravitational field of a body outside of its matter distribution in its own 
local coordinate system, and improved translational and rotational equations of motion for 
a system of A^-extended rotating gravitationally interacting bodies of arbitrary shape and 
composition. This framework, in our opinion, shifts the 'pseudo Newtonian view' of the 
classical PPN-framework towards a more relativistic one. Note, that this point already was 
the motivation for the DSX-formalism. All main results from the classical PPN-framework 
as mentioned in the introduction (Section |I[) are recovered, but now with the help of an im- 
proved and more consistent framework. Not only the usual PPN Einstein-Infeld-Hoffmann 
equations of motion for a system of mass monopoles, the spin-orbit and spin-spin interaction 
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terms have been derived in the new (3, 7-framework, but also all the various terms (depend- 
ing upon ?7, (7 — 1) or [j3 — 1)) related with a violation of the Strong Equivalence Principle 
are recovered. 

It is the hope that this version of a PPN-formalism might eventually form the basis of 
a new test-theory for experimental gravity. To this end more work has to be done mainly 
concerning: a discussion of orders of magnitude of the various terms that appear in expan- 
sions in many places, the construction of models for the astronomical bodies, e.g., for the 
determination of the time dependence of the various mass and spin multipole moments of 
the bodies and, last not least, the relation of formalism with observations. This last point 
does not only concern the theoretical formulation of the various relevant measuring tech- 
niques (VLBI, SLR, LLR, GPS, astrometry, etc.) by employing the new PPN- framework, 
but the usefulness of certain theoretical concepts introduced here, e.g., the post-Newtonian 
spin of a body, has to be discussed. It is obvious that one wants to keep concepts that 
have been extremely useful in Newtonian classical mechanics for reasons of simplification 
also in a relativistic framework. However, the fact that many of these concepts do not even 
have a meaning in GRT (without resorting to some approximations) for real astronomical 
situations sheds some light on the difficulties that one faces here. On the other hand it is 
also obvious that the relevant measuring techniques at least in certain fields of application 
such as the problem of Earth's rotation have reached an accuracy level where a wealth of 
physical effects enters that simply cannot be modeled. In that case it has become a practice 
to start with a theoretical 'approximate model' and to measure 'offsets' from the theoretical 
predictions. What parts of the new PPN-framework presented here will eventually play a 
role in such an 'approximate model' (e.g., for the description of Earth's rotation) in order 
to reduce residuals is not so clear at present. 
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